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^ ! In this paper, we use various ansatzes with undetermined functions and the technique of 
• 1— I ' 

c/3 ■ moving frame to find solutions with parameter functions modulo the Lie point symmetries for 
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i-C ! the classical non-steady boundary layer problems. These parameter functions enable one to find 
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the solutions of some related practical models and boundary value problems. 
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22 ■ 1 Introduction 



In 1904, Prandtl observed that in the flow of slightly viscous fluid past bodies, the frictional 

effects are confined to a thin layer of fluid adjacent to the surface of the body. Moreover, 

- he showed that the motion of a small amount of fluid in this boundary layer decides 

H ' such important matters as the frictional drag, heat transfer, and transfer of momentum 
cz ' 

between the body and the fluid. The two-dimensional classical non-steady boundary layer 
equations 

Ut + UUx+VUy+Px=Uyy, (l.l) 

Py = 0, U:r + Vy = (1.2) 

are used to describe the motion of a flat plate with the incident flow parallel to the plate 
and directed to along the x-axis in the Cartesian coordinates {x,y), where u and v are 
the longitudinal and the transverse components of the velocity, and p is the pressure. Lie 
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point symmetries of the above equations were obtained by Ovsiannikov [O]. The three- 
dimensional classical non-steady boundary layer equations are: 

Ut + UUx + VUy + WUz^~Px + T^Uyy, (1.3) 

Wf + UWx + VWy + WWz — Pz + ^Wyy^ (l-^) 

Py = 0, U:^ + Vy + Wz = 0, (1.5) 

where p is the density constant and v is the coefficient constant of the kinematic viscosity. 
The Lie point symmetries of the equations (1.3)-(1.5) were obtained by Garaev [G] and by 
Vereschagina [V]. Lie group method is one of most important ways of solving nonlinear 
partial differential equations. However, the method only enables one to obtain certain 
special solutions. It is desirable to find more effective ways of solving nonlinear partial 
differential equations. 

The aim of this paper is to develop more powerful methods than the Lie group method 
of solving the above classical non-steady boundary layer problems. Based our earlier work 
[X] on function-parameter exact solutions of the equation of nonstationary transonic gas 
flow discovered by Lin, Reisner and Tsien [LRT], and its three-dimensional generalization, 
we find ansatzes with various undetermined functions and obtained solutions with param- 
eter functions modulo the Lie point symmetries for the above boundary layer equations. 
In the three-dimensional case, we have also used the technique of moving frame. Our 
results can be applied to certain related practical models and boundary value problems 
by specifying these parameter functions. Below wc give a more detailed introduction. 

Lin, Reisner and Tsien [LRT] found the equation 

2Utx + UxUxX - Uyy = (1.6) 

for two-dimensional non-steady motion of a slender body in a compressible fiuid, which 
was later called the "equation of nonstationary transonic gas flows" (cf. [M]). Using 
certain finitc-dimcnsional stable range of the nonlinear term UxUxx, we obtain in [X] the 
following solution of the above equation: 

where /3 is an arbitrary function of t. Since the above solution blows up on the moving 
line + /3 = 0, it reflects partial phenomena of gust. Moreover, we have found a 
family of solutions of the three-dimensional generalization of the equation (1.6) blowing 
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up on a rotating and translating plane y cos a{t) + z sin a{t) — f{t), which reflect partial 
phenomena of turbulence. Our work [X] suggested some possible approaches to more 
general nonlinear partial differential equations in fluid dynamics, such as the classical 
non-steady boundary layer problems we will deal with in this paper. 

By the divergence-free condition in (1.2), we write u and v in the potential form: 

u^^y, (1.8) 

where ^ is a function of t, x, y. It is known that the solution space of the equations (1.1) 
and (1.2) is invariant under the following transformation T (cf. [O]): 

T{u)^u{t,x,y + a), T{v) ^ v{t,x,y + a) - at - a^u, T{p)^p (1.9) 

for any function a of t, x. Modulo this transformation, our main method of solving the 
equations (1.1) and (1.2) is to assume 

i = fit, x) + K{t, x)y + g{t, x)H{t, y), (1.10) 

where /, k, g are functions to be determined and if is a more sophisticated function that 
is either first given or is to be solved under certain conditions. The term K,{t, x)y is 
used to make the function H more flexible. With the ansatz (1.10), the equation (1.1) 
is reduced to a system of several nonlinear partial differential equations in t, x. While 
by Lie group method, one usually reduces it to a single nonlinear partial differential 
equation in t, x. This in a way shows that our method is more powerful than the Lie 
group method. Imposing some conditions on the above undetermined functions, we are 
able to solve the system of partial differential equations. In this way, we obtain three 
families of function-parameter exact solutions of the two-dimensional classical non-steady 
boundary layer problem (1.1) and (1.2). For instance, we have the following solutions (see 
Theorem 3.2): 

^ = ^ + ^» cos(V^Z/), (1.11) 



2/3' ^ 

^ = - ) MVP' y), (1-12) 



_ {f3"f - 2/3' 13"' 2_±-20 
^ 4(/302 2/3''' 



X 



(1.13) 



where /5 is an arbitrary increasing function of t and 3" is any one-variable function. Ap- 
plying the symmetry transformations of the equations (1.1) and (1.2) (see next section) 
such as the one in (1.9), we can get solutions with more parameter functions. 
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The main contents of this paper are devoted to solving the three-dimensional clas- 
sical non-steady boundary layer equations (1.3)-(1.5). Similarly, we write u,v,w in the 
potential form: 

u = ^y, w = r)y, v=-{^^ + r],) (1.4) 

for some functions ^ and r] in t, x, y, z. The equations (1.3)-(1.6) have the similar symmetry 
transformation as that in (1.9). Modulo this transformation, we assume 

C = /(t, X, z) + nit, X, z)y + g{t, x, z)H{t, a{t, x, z)y), (1.14) 

and 

rj = r{t, x, z)y + ({t, x, z)^{t, a{t, x, z)y), (1.15) 

where /, g, k, a, r, ( are three- variable functions to be determined, and H{t, zu) and $(t, zu) 
are more sophisticated two- variable functions that are either first given or are to be solved 
under certain conditions. Again the terms K{t,x,z)y and T(t,x,z)y are used to make H 
and $ more flexible. With the ansatz (1.14) and (1.15), the equations (1.3) and (1.4) are 
reduced to a system of several more complicated nonlinear partial differential equations 
in t, X, z. This again in a way shows that our method is much more powerful than the 
Lie group method. One of the main technique of solving the system of nonlinear partial 
differential equations in t, x, z is to use time-dependent orthogonal transformations (mov- 
ing frames). Imposing some technical conditions on the above undetermined functions, 
we are able to solve the system of partial differential equations. In this way, we obtain 
ten families of function-parameter exact solutions of the three-dimensional classical non- 
steady boundary layer problem (1.3)-(1.5). For instance, we have the following solutions: 
(1) for any function i of z, any function a of t and arbitrary functions '<s,(p in w = a — l: 

dz 



ii 

a if 

— TTT^ F= 

2q;' J^. 



exp (^2u^'^ J 



at y 




cos I 1 , (1.16) 



(6.17) 



P--I[-^(^^'^W^- (6-18) 
(see Theorem 6.1); (2) in terms of the moving frame 

w = xcosj + zsinj, w = zcosj — xsinj, (1-19) 

l3^-,' + 2;/(.r cos - - 23- ' sin - - 67') , 

u — : — 7- r (tu -|- 2ti7tan7) 

4i/(/3-6i/sin7) ^ 

— ^^/^7' 7 + 7' 7^ -07 — {Gux + Puj)y~'^, (1-20) 



/3^y + 2v{P' COS 7 - 2Pi sin 7 - 6y ) /^Y 

= 77, — 7 — : ^ y - sec 7 - 6z/|/ , 1.21 

4i/ cos 7 (p — Di/ sm 7j Di/ 

/?2y + 2ui(3' COS 7 - 2l3i sin 7 - 67') , 

w = — ^ — - — : r (rotan7-2ro) 

4z/(p — 6i/sin7) 

+^/97'ti7 — 7'rosec7 — {Qvx + (3w)y~'^ t&n.^ (1-22) 

and p is more complicated and given in (7.26), where (3 and 7 are arbitrary functions of 
t (see Theorem 7.1). Applying the symmetry transformations of the equations (1.3)-(1.5) 
(see next section) such as that similar to the one in (1.9), we can get solutions with more 
parameter functions. For instance, the above second solution would yield to a solution 
singular on any moving hypersurface, which may be used to study turbulence. In addition, 
we present another two anzatzes and obtain another two families of function-parameter 
exact solutions of the three-dimensional problem that are rational functions in y. 

Of course, one can view the two-dimensional problem (1.1) and (1.2) as a special 
case of the three-dimensional problem (1.3)-(1.5). Our results on the two-dimensional 
problem in Section 3 can be directly applied to the three dimensional problem by assuming 
that ^ is independent of z and r] is independent of x. Applying constant orthogonal 
transformations (see Section 2), we can obtain more sophisticated solutions of the three- 
dimensional problem. So obtained solutions of the three-dimensional problem will not be 
given. Thus not all the solutions of the two-dimensional problem in Section 3 are the 
specializations of the solutions of the three-dimensional problem in later sections. Our 
another purpose of presenting the solutions of the two-dimensional problem in Section 3 
is to let the reader first understand our techniques in simpler version and then to better 
understand much more sophisticated techniques of solving the three-dimensional problem 
in later sections. Indeed, due to the technicality, the author himself keeps referring the 
results and arguments on the two-dimensional problem in Section 3 when he solves the 
three-dimensional problem in later sections. 

For convenience, we always assume that all the involved partial derivatives of related 
functions always exist and we can change orders of taking partial derivatives. As we all 
know that in general, it is impossible to find all the solutions of nonlinear partial differ- 
ential equations both analytically and algebraically. In our arguments throughout this 
paper, we always search for reasonable sufficient conditions of obtaining exact solutions. 
For instance, we treat nonzero functions as "nonzero constants" for this purpose because 
our approaches in this paper are completely algebraic. Of course, one can used our meth- 
ods in this paper to get more solutions, in particular, by considering the support and 
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discontinuity of the related functions. 

The paper is organized as follows. In Section 2, we present the known Lie point 
symmetries of the equations (1.1) and (1.2) due to Ovsiannikov [O], and the Lie point 
symmetries of the equations (1.3)- (1.5) due to Garaev [G] and by Vereschagina [V]. In 
later sections, we always try to exclude the obvious ingredients in our solutions induced by 
these symmetry transformations. In other words, we search for the solutions modulo these 
transformations. We solve the two-dimensional classical non-steady boundary layer equa- 
tions (1.1) and (1.2) in Section 3. Starting from Section 4, we solve the three-dimensional 
classical non-steady boundary layer equations (1.3)-(1.5) according to the types of the 
involved functions with respect to the variable y. First in Section 4, we find solutions 
of the equations (1.3)-(1.5) containing a single exponential function in y. In Section 5, 
we look for solutions of the equations (1.3)- (1.5) with distinct exponential functions in 
y. In Section 6, we obtain solutions of the equations (1.3)-(1.5) with trigonometric and 
hyperbolic-type functions in y. Finally in Section 7, wc present solutions the equations 
(1.3)-(1.5) which are rational functions with respect to the variable y. 

2 Lie Point Symmetries 

In this section, we will present all the known Lie point symmetries of concerned boundary 
layer problems. 

Denote by T(i){xi, ...,Xn) the space of functions in xi, whose all partial deriva- 

tives up ith order exists. Ovsiannikov [O] found the following symmetric group G2 gen- 
erated by the transformations 

{^la, T2b, Tab, T4^, T^f), TG^\a,beR,b^O;ae T(2)(t)\l3 e T{o)(t)]a e T(i)(t, x)} (2.1) 

with 

Tia{u) = u{t + at,x,y), Tiaiv) ^ v{t + at,x,y), Tia{p) = p{t + at,x); (2.2) 

T2b(u) = u(b\ b^x, by), T2b(v) = b-^v(bH, b^x, by), T2b(p) = p(bH, b^x); (2.3) 

T-ib{u) = bu{t, bx, by), T^b{v) = v{t, bx, by), T^b{p) = b^pit, bx); (2.4) 

T4a(u) ^u(t,x + a,y) —a', T4a(v) ^ v(t,x + a,y), T4a(p) ^ p(t,x + a) + a"x; (2.5) 

T5p{u) = u, Tr,p{v) = V, T^oip) =P + (2.6) 

Teaiu) ^u{t,x,y + a), Ta^iv) ^ v{t,x,y + a) - at - a^^u, T^^ip) ^ p. (2.7) 
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The elements oi G2 transform the solutions of the two-dimensional classical non-steady 
boundary layer equations (1.1) and (1.2) into solutions. For instance, the transformation 
T^aT^ transforms a solution with zero pressure into a solution with the pressure p — 
a''x + 7 and the transformation Tqo^ may change a solution to the one with an additional 
two-variable parameter function. In this way, we can always transform our solutions to 
the ones with nonzero pressure and at least one two-variable parameter function. For 
simplicity, we only write our solutions of (1.1) and (1.2) as simple as possible modulo the 
group G2- We will do the same thing for the solutions of the three dimensional problems 
(1.3)-(1.5). 

The symmetry group G3 of the equations (1.3)-(1.5) are generated by the transforma- 
tions 

{^lo, T'2a, T^b, Tib, T^a, Tea, Trp, T^a I a, 6 e M, 6 7^ 0; 

a e J^(2)(t);(3 e J^(o)(0; e ^mit, x, z)} (2.8) 

with 

TiaW ^i/jit + a^x^z^y), i/j ^ u,v,w,p; (2.9) 

T2aiu) — u{t,xcosa + zsma,—xsma + zcosa,y)cosa 

+w{t,xcosa + zsina, —a; sin a -|- zcosa,y) sin a, (2-10) 

T2a{'w) — — u{t,x cosa + zsina, —xsina + cosa,y) sina 

+w{t,x cos a + 2; sin a, — xsin a + zcosa, y) cosa, (2-11) 

T2a{i^) — i^{t.iXCOsa-\- zsina,— xsva.a-\- zcosa,y), i^ — '^iP', (2-12) 

Tzb{u) = bu{t, bx, y, bz), T3b{v) = v{t, bx, y, bz), (2-13) 

Tsbiw) = bw{t, bx, y, bz), TM = b''p{t, bx, y, bz); (2.14) 

Tib{u) =b~'^u{bH,x,by,z), T^hiv) = b~^v{bH,x,by, z), (2.15) 

T^biw) = b-^w{bH, X, by, z), T^bip) = b-^p{bH, x, by, z); (2.16) 

T5a{u) = u{t,x + a,y,z)-a', Tr^^{v) = v{t,x + a,y,z), (2.17) 

T^a{w) = w{t,x^ a,y,z), T^^ip) = p{t,x + a, z) + pa"x; (2.18) 

Tfiaiw) ^ w{t,x,y,z + a) - a', Ta^iv) ^ v{t,x,y, z + a), (2.19) 

Teaiu) = u{t,x,y,z + a), Ta^ip) = p{t,x, z + a) + pa" z; (2.20) 



Ty0{u)^u, T,p{v)^v, Tjf,{w)^w, Trp{p)^p + P; (2.21) 

Tg^iu) = u{t,x,y + a,z), Tg^iv) = v{t,x,y + a) - at - a^u - a^v, (2.22) 

Tsa(u) ^w(t,x,y + a, z), T^aij)) = p. (2.23) 

3 Solutions for the 2-D Problem 

In this section, we will find various function-parameter exact solutions for the two- 
dimensional classical non-steady boundary layer equations (1.1) and (1.2). 
Solving the second equation in (1.2), we get 

u^^y, -^^ (3.1) 
for some function ^ in t,x,y. Now the equation (1.1) becomes 

^yt ~l~ ^y^xy ^x^yy ~l~ ^yyy ("^-2) 

Assume that ^ is of the form: 

C = f{t, x) + n{t, x)y + g{t, x)H{t, y), (3.3) 
where K{t,x), g{t,x), f{t,x) and H{t,y) are two- variable functions. Then 

^y = K + gHy, ix = fx + K'xV + QxH, iyt = Kt + 9tHy + gHyt, (3.4) 

^yx ~l~ 9xHy, ^yy gHyy, ^yyy d^yyy (3-5) 

In particular, the nonlinear term in (3.2) becomes 

^y^xy — ^x^yy = (^ + gHy){Kx + g-x Hy) - {fx + Kxy + gxH)gHyy 

= KK^ + {Kg,; + K,;g)Hy - {f^ + K^y)gHyy + ggx[Hy - HHyy]. (3.6) 

In order to linearize (3.6) with respect to H, we divided it into the following cases. 
First we consider H = e^^ . Then 

^yt = «t + {iQt + 7V)e'^" + TT^ye^^. (3.7) 

So (3.2) becomes 

nt + ni^x + b'g + i{gt + I'gy + i^Qx + Hxg - ifx - ii^xgy - i'^9)]e''" + = o, (3.8) 

which is implied by the following system of partial differential equations: 

Kt + KKx +Px = 0, 7' - 'JKx = 0, (3.9) 



ig + iQt + K-iGx + iK-xQ - i^g - I'^gfx = o. (3.10) 

So 

«;=^X + 7o (3.11) 

7 

for some function 70 of Moreover, (3.10) implies 

/a; = ^ H ^ ^ 7- 3.12) 

7 75' 19 



So 



e = / + ^xy + 7oy + ^e^^. (3.13) 
7 



Hence, 

^-t/-'-^^ (3.14) 

modulo some T^p in (2.6). 

Next we consider the case f — Qx — 0- Replacing H by gH, we can assume g — 1- 
Now (3.2) becomes the following linear partial differential equation: 

Kt + HKx + Hyf + KxHy Kxl/Hyy + ^j; = H yyy . (3.15) 

To simplify the problem, we assume 

7" 

Kt + K/ta; + = 0, Kx^— (3.16) 



for some nonzero function 7' of Then 

~ 27 

for some function 70 of t and 



7" 

K = —X + 7o (3-17) 



27r-(y? ^. , 7oy^+7^y ^ .3,g. 

modulo the transformation in (2.6). Now (3.15) becomes 

y 

Hyt + ^(-^2/ ~ ^^w) = Hyyy. (3.19) 

Write 

Hy = -^H{t,^y) (3.20) 
for some two variable function H{t, w) with = ])■ Then (3.19) turns out to be 

iit = I'H^^. (3.21) 
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One can use Fourier expansion to solve the above equation with given initial condition. In 
particular, discontinuous solutions can be found in this way. For the algebraic neatness, 
we just want to find globally analytic solution with respect to the spacial variables x and 
y. So we take the following solution of (3.21): 

m 

H^Yl (iie("?-^?)^+"'^ sin(6,(2a,7 + zu) + a), (3.22) 

i=l 

where Oj, bi, q, di are real numbers. Hence, we have 

^ m 

= -^5]^^^^"^"'''^'"'"'^'"^"(^^(2a,7+ y) + cO. (3.23) 

By (3.1), we have the following result, where the first solution follows from a direct 
observation the equations in (1.1) and (1.2). 

Theorem 3.1. We have the following solutions of the two-dimensional classical non- 
steady boundary layer problem (1-1) and (1-2): (1) 

u^y, v^-Kj:, P^k, (3.24) 

where n is an arbitrary function oft,x; (2) 

= — + 7o + 7^e^", (3.25) 
7 

V = 7 J g^e^y. (3.26) 

7 19 19 1 

and p is given in ( 3. 14 ), where 7, 70 are any nonzero functions of t and g is any nonzero 
function oft,x; (3) 



7 X 

u 



^ m 

^ + 70 + ^ V rf.e('^' -''?)^+»^^^ sin(6,(2a,7 + v7 ?/) + q), (3.27) 
27 V7 ~( 

V = -0 (3.28) 

and p is given in (3.18), where 7,70 are any functions of t, and ai,bi,Ci,di are real 
constants. 

Next wc consider the case / = dy{Hy — HHyy) = 0. Again let 7 be a function of t. 
For a e M, we denote 

^0 = ^1 = sin(7y), (3.29) 



^0 = -1^ , ^1 = cos(7|/). (3.30) 
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Given r e {0, 1}. Assume 



Then 



(3.31) 
(3.32) 



Hi - HHyy = (5o,ra + KrW- 
By (3.6), 

iyiyx - ixiyy = K-K-x + {^-Qx + K'x9)l^r " (-l)''7^(/x + '^xy)g'&r + (^Cr^ + h,rW99x- (3.33) 

Since 

= + (5*7 + 9l'¥r - {-lyiigy^r. iyyy = {-iTl^, (3-34) 



Kt + + Px + +b'g + i{9t + i^gx + Ka^s- - {-iy'j'^g)]dr 

-{-lY'fi'y' + 'yKx)yg^r + (5o,ra + Si,rh'^ggx = 0, 
which is imphed by the following equations: 

Kt + + {So,ra + Si^rWgg^ +Px^^, li^x + 7' = 0, 

gi + {gt + ngx + /^xs- - (-i)''7^5')7 = o. 

For convenience of solving the above equations, we assume 

A37 ^ y /3" 



So we take 



K 



7 



g 



7 2/3' 

e(-i)''/3 



2/3' ' ^ /?' 
for some one-variable function Q'. Moreover, 



X 



^ 4(/3')2 ^ 2/3' 

modulo the transformation in (2.6) and 



X 

7w 



_ (3"xy e^-^yi^ 
2/3' ^ /3' 



a; 



'&r. 



(3.35) 

(3.36) 
(3.37) 

(3.38) 
(3.39) 

(3.40) 
(3.41) 



According to (3.1), we have: 



Theorem 3.2. Let j3 he any nonzero increasing function oft and let $5 he any one- 
variahle function. We have the following solutions of the two-dimensional classical non- 
steady boundary layer problem (1.1) and (1.2): (1) 



u 



X 



p"x 
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(3.42) 



y = 2^^' ( ) ^e^y _ ae-^y) (3.43) 



2/?' 2y(^ vv^ 

and p is given in (3.40) with r — 0, where a e R; (2) 

(3"x e-^ 



u 



,;T^'^^{-^)^^<^'y)^ (3.44) 



v^-^- -4^'^' ( ] sm(Jp' y) (3.45) 



anc?p is given in (3.40) with r — 1. 

Finally, we look for ^ of the following form: 

e = K{t, x)y + g{t, x) + h{t, x)y-\ (3.46) 
where g^, /i are functions of x with 0. Then 

Cj/ = « - ^J/t = «t - ^ty~^ Cyx = - ^a;y~^ (3.47) 

= «a;y + + h^y~^, iyy = ^2/2/2/ = -Qhy~^. (3.48) 

So (3.2) becomes 

Kt - hty'"^ + (k - hy~'^){K^ - h^y'"^) - 2hy~^{Kxy + gx + h^^y"^) +Px^ -^hy"^, (3.49) 
equivalently, 

Kt + +Px- {ht + kHx + ?>Kxh)y~'^ - 2hgxy~^ - hh^y'^ = -Qhy'*. (3.50) 
Moreover, (3.50) is implied by the following equations: 

Kt + KKx + Pa; = 0, ht + kHx + 2>Kxh — 0, (3.51) 

^ = 0, hx^Q. (3.52) 

So /i = 6a; modulo T^a in (2.5). 

According to the second equation in (3.51), 

6k + l^xKx = =^ K = 7x~^/^ (3.53) 

for some function 7 of Thus 

^ = 7x-^/3^ + Qxy-^. (3.54) 
12 



Moreover, Kt — ^'x Thus 

by modulo T^p in (2.6) by the first equation in (3.51). Therefore, (3.1) implies: 

Theorem 3.3. Let 7 he any function of t. We have the following solutions of the 
two-dimensional classical non-steady boundary layer problem (1.1) and (1.2): 

u^^x-'/'-6xy-', v = lx-'/'-6y-\ p = ^x'/' - ^x-'/\ (3.56) 



4 3-D Uniform Exponential Approach 

In this section, we will find certain function-parameter exact solutions with single ex- 
ponential function in y for the three-dimensional non-steady boundary layer equations 
(1.3)-(1.5). 

Solving the second equation in (1.5), we get 

U^^y, W^riy, -(^a: + Vz) (4.1) 

for some functions ^ and rj in t, x, y, z. Now the equation (1.3) and (1.4) become 

Cyt + CyCyx — {Cx + Vz)Cyy + VyCyz + -Px = ^Cyyyj 

(4.2) 

Vyt + CyVyx - {Cx + Vz)Vyy + VyVyz + -Pz = ^Vyyy. (4.3) 
First we assume that ^ and 77 are of the form: 

i = f{t, X, z) + n{t, X, z)y + g{t, x, z)H{t, a{t, x, z)y), (4.4) 

and 

T] = T{t, X, z)y + C{t, X, z)^{t, a(t, x, z)y), (4.5) 

where /, g, k, a, r, ( arc three- variable functions, and H{t, w) and $(t, w) are two variable 
functions, where zu = ay. Then 

^y = n + agH^, r]y^T + aC^^, (4.6) 

^x^ fx + K,xy + QxH + a^gyH^, rj^ = r^y + (z^ + cr^C^^, (4.7) 
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^yx ^ K'x + {(yg)xH^ + CrCTxgyH^-cj, Vyx ^ + {ctOx^vj + CTCTxgy^vju,, (4.8) 

^yz = >iz + {(Tg)zH^ + aa^gyH^^, rjy^ = + {crC)z^^ + cra^^gy^^^, (4.9) 

^yy ^ fl'-^roro; Cyyy ^ fl'-^rororoj T]yy C^roro; Vyyy ^ C^zazu-uc- (4-10) 

Thus the nonhnear term in (4.2) becomes 

^y^yx - {^x + Vz)^yy + Vy^yz = (« + (7gH^){n^ + {ag)^H^ + aa^gyH^^) 

+ {^^x + Tz)y + gxH + (T^gyH^ + + cr;j(^|/$^) 
+ (r + (TC$^)(fi;^ + {ag);,H^ + aa^gyH^^) 
= KKj; + TK^ + [{crgK)^ + {ag)^T]H^ + (tC,k^^vj + (7g{[K,(Ta: + Tcr^ - a{Ka: + T^)]?/ 

-(yfx}H^v. + (Jg{ag)^Hl - a'^gg^HH^^ + aC{ag),H^^^ - a'^g^H^^^. (4.11) 

Symmetrically, the nonlinear term in (4.3) becomes 

CyVyx ~ {Cx + 'nz)Vyy + VyVyz 
= KT^ + TT^ + [((jCt)^ + {aQ^K]^^ + crgr^H^ + crCHKa^ + ra^ - a{Ka: + T^)]y 

-crfx}^^^ + crC{(TC)z^l - cr'^CCz^^^z^ + crg{aC)xH^^^ - a'^g^CH^^^- (4-12) 

In order to linearize (4.11) and (4.12) with respect to H and we divide it into ten 
cases. 

Case 1. H = ^ = e^. 

In this case, 

iyiyx - {^x + 'nz)iyy + Vy^^y^ = H^x + TH^ + {(TgH^ + uC^z + {'7g)xK' + ((7fi')^T 

+(Tg{[n(T:, + rcr^ - cr(fi;a; + Tz)]y - af^jje"^ + a^j^g"^ + C,{(jg)z - agQe^"^ . (4.13) 
Moreover, 

^yt^ Kt + (((Tg)t + (T(Ttgy)e'^, ^yyy ^ (y^ ge"^ . (4.14) 

Thus (4.2) becomes 

Kt + KK^ + TK^ + {((T5f)t + CrgK^ + (tCk^ + {(Tg)xK- + ((T5f)2T + (T5f{[(Tt + Ka^ + T(T^ 

-ct(/^^ + T,)]y - - i^a^}}e^ + a[a^g^ + ((a^). - agCz]e^^ + -Pa. = 0. (4.15) 

P 

By symmetry, (4.2) and (4.3) are imphed by the following system of partial differential 
equations: 

Kt + KKx + TKz + -p^ = 0, rt + KT^+TTz + -Pz = 0, (4.16) 

P P 
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at + Ka^ + ra;, - a(K,x + rz) ^0, (4-17) 

{ag)t + agn^ + aCn^ + {ag).^K + {ag)^T - a^g{fx + = 0, (4.18) 

{aQt + (TQT^ + (tCtz + ((7C)x« + {(jOzT - t'^CUx + 1^(7) = 0, (4.19) 

(^x/ + C(f^^7). - t^^C. = 0, + 9H)x-(ygxC^^- (4.20) 
Observe that (4.20) yields 

O-a; + (T^- = (T ( - ) , (Tz + (^x^ = ( ^] ■ (4.21) 



Thus we take 

j). (^) (c)^ ^ (^) (^)z i.^yx 
Modulo T^a^Tea in (2.17)-(2.20), we take the solutions 



(4.22) 



C z 

- = tan7, -, (4.23) 
9 X 

where 7 is a function of t. In the first case, 

CTj; + CT^ tan 7 = 0, (4.24) 

and we take 

^ _ gO(2COS7-=.sin7)gg^^ (4.25) 

for later convenience. In the second case, 

xax + zaz = cr, (4.26) 

and we take 

(4.27) 

where 71 is a function of t and 3" is a one-variable function. Furthermore, ag x (4.19) — 
cr( X (4.18) gives 

(aC)tcrfi' - aCiag)t + {orgyr^ - {aCfn^ + (aC) (crfl-) (r^ - k^) 

+ [(<7C)x<75 - (^C((^5)J/« + [H)z'y9 - c^C(<75)Jt = 0. (4.28) 
Multiplying (4.28) by l/^agf, we obtain 

Since Pa;2 by our convention, (4.16) implies 

(«2 - Ta;)t + (/^(k^ - Tx))x + ('r(«;2 - ra;))^ = 0. (4.30) 
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Assume CI 9 — tan 7. Denote the moving frame 

= XCOS7 + 2;sin7, = 2;cos7 — xsin7. (4-31) 

Then 

— cos ^ dx + sin 7 dz, — — sin ^ dx + cos 7 dz- (4.32) 

Observe that (4.29) becomes 

' "Tx ~ K'z tan^ 7 + (Tz — Kj;) tan 7 = 0, (4.33) 



cos-^ 7 



equivalently, 

7' + (cos 7 + sin 7 92)(tcos7 — k sin 7) = 0^7' + d^{Tcosj — fi;sin7) = 0, (4.34) 
and so we take 

T COS 7 — K sin 7 = — 7'ro =^ T — K tan 7 — 7'ro sec 7 (4.35) 
for the convenience of solving the problem. Note 

dt{w) — 7'g7, dt{'dj) = —7' 07. (4.36) 
Substituting (4.25) and (4.35) into (4.17), we obtain 

27' sin 7 - 27' cos 7 w^' - d^{n) = 0, (4.37) 
where $5 is a function of w. Thus we take 

K = 27'tz7 sin 7 — 7'ot^Q'' cos 7 + £(t, (4.38) 
for some two- variable function e{t,w). Moreover, 

COS 2^ 

T = —^'w ^'vj'^^' sin 7 + £ tan 7. (4.39) 

cos 7 

Note 

«2 - = -7'ro^$j" + sec 7 + 7'. (4.40) 

Moreover, 

(«. - r,)t = -7"^'$>" - 2{ifww^" + (7')'^'^'" 

+ sec 7 (e^t + 7' tan 75^- 7'tz7£^^) + 7", (4.41) 
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- T.)]. + [t(«, - = {if [a'-s'^" + ^"'\ - 3(7')' tan 7 
X ^"w^ - 7'[sec 7(2£$>' + e^) + 27' S]^^ + 7' tan 7 (sec 7 £^ + 7')- (4.42) 

Since the subcase 7 e R is included in Case 3 modulo the transformation in (2.10)-(2.12), 
we assume 7' ^ 0. Thus (4.30) is implied by the following system of partial differential 
equations: 

+ 2$J'$J" = 0, (7" + 3(7')' tan7)»" = 0, (4.43) 
7'' + (7')^tan7 + 27'£^sec7 tan7 + £^tsec7 = 0, (4.44) 
[(£9' + e^) sec 7 + y^^^l^ = 0. (4.45) 

By (4.44), we take 

e — iiw') cos^ 7 — iw cos 7 (4.46) 
for some function l of w. Moreover, (4.45) is imphed by 

(^J'(n7)t(tJ7) + t'(tJ7))^ = 0. (4.47) 

First, 

7 is arbitrary and t = ci + d\e~"''^ if = aw + 6, (4.48) 
where a,b,Ci, di e R. Assume 7^ 0, we have 

^'"{w) + 2»'(OT)a"(OT) = 0, 7" + 3(7')^ tan 7 = 0. (4.49) 

Thus 

S> = \n(aw + 6), lnsin(aro + 6), ln(6e'^^ + ^le""^) (4.50) 
with a, 6, 61 e M and 7 is implicitly given by 

= Ce'^t-2sec7 tan7^ c, (i G R. (4.51) 

1 — sin 7 

Furthermore, 

L = ci(aro + 6) + , if = ln(aro + 6), (4.52) 

t = Cl cot(aro + 6) + (ii csc(a-^ + 6) if = lnsin(ati7 + 6), (4.53) 

i^ci-^. ^i— ^ + — — ^ ^ if $J = ln(6e"^ + 6ie-»^). (4.54) 

According to (4.35), (4.38) and (4.39), 

— Kt + K.d^{K,) sec 7 — iwKz sec 7 

= 2(7'' sin 7 + (7')^ cos 7)ot + 2(7')^ti7 sin 7 + — ^ws^ 
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+((7' )^ sin 7 - 7" cos -f)w'^Q' - 2(-f'yw'njQ' cos 7 + (7')^^^$^" cos 7 
+ (27'rosin7 — 7'-ji7^S>' cos 7 + (27' tan 7 — 2j'zu'^') 
—'y'w{2'y' sin^ 7 sec 7 — 7'(2tz/Q'' sin 7 + tz7^Q='' cos 7) + e^) 
= 2(7')^^^(9" + (9')^)cos7- (3(7')^sin7 + 7"cos7)^z72$>' + £t + 27'(£tan7 

+7'rosin7) + 2(7'' sin 7 + (7')^ sec 7 — 7' (7'^$^' cos 7 + £^ + ^'s))'^, (4.55) 

— Tt + Kd^{T) sec 7 — j'wTz sec 7 

/x9 . / ,x 9 COS 27 sin 7 ,/COs27\ ^ , ,,nCOs27 ^ ^, ,,9^ . 

= 4 y f sin 7 - 7'^ ^ - 7 ^ ^ - 7' - 2 y f OTG7^'sin7 

\ cos^7 cos 7 / cos 7 

— (7''sin7 + (7')^cos7)^^$J' + (7')^^^$>''sin7 + 7'sec^ 7 £ + (st — 7'^ e^) tan 7 

(cos 27 \ 
7^ h 2tzjQ='tan7 J + 7'^^' cos 27 tan 7 
cos^ 7 / 

+^'w(2Q' sin^ 7 + wQ'' sin 7 cos 7) — £^ sin 7] sec 7 
= 2(y)^^^(^>"+ (^>')^)sin7- (3(y)^sin7 + y'cos7)ro^$5'tan7 
+27'(7'sec7 — (7'w53'cos7 + + $5'e))tz7tan7 
^"^^ (y'^ + (y + (£t + 2e tan 7) tan 7. (4.56) 



cos 7 
Note 

C in Cases (4.48) and (4.54), 

^" + (cj')2 = 5 ^ J in Case (4.52), (4.57) 
[ — in Case (4.53), 

{aci in Cases (4.48) and (4.54), 
2aci in Cases (4.52), (4.58) 
—aci in Case (4.53), 

and 

£t + 2£ tan 7 = (7')^^z7cos7 — ((7')^ sin 7 + 7'' cos 7)117 (4.59) 
by (4.46). Moreover, 

2^X7 sin 7 — ■;S7COS7 — 3a; sin 7 cos 7 + 2;(3sin^7 — 1), (4.60) 

-07 — ■zi7sin7 = x(3sin 7 — 1) + 2;(1 — 3cos 7) tan 7, (4-61) 



cos 7 

cos 27 



— sin 7 tan 7 = — cos 7. (4.62) 



cos 7 

By (4.16), (4.43) and (4.56)-(4.63), we have: 

,3(y)Han7 + y^3^, {-f')^[w^ - {bw^ + ,.2,^ 
p = p{^^ — ^-^ —uj'^^'+^:-^^ )- — ^ + y'ca^(c- 7 sec 7) sec 7 

,/ Q • 2 N 3a;2sin27 ^2(3^^32 ^^^g^^^ x 
+7 ( a;2;(l -3sin 7) + 1} (4.63) 
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modulo the transformation in (1.26). 

Modulo the transformation given in (2.22) and (2.23), we take 

g ^ ^ cos-i e-^^^^ (4.64) 

by (4.25). According to (4.18), (4.23) and (4.38), we have 

fx = 27'e~^(sin7 — w'^' cos'~)) — ve^ secy. (4.65) 

Furthermore, 

^ — f + (27'rosin7 — ^'w'^Q' cos^ + e)y + cos 7 exp(e^ysec7 — 3^), (4.66) 

f COS 2^ \ 

77 = ( £tan7 — 7' w — 7'ro^Q^'sin7 1 y + sin 7 exp(sec7 e^y — (4.67) 

\ cos 7 / 

By (4.1), we have the following theorem. 

Theorem 4.1. In terms of the notations in (4-31), we have the following solutions 

of the three-dimensional classical non-steady boundary layer equations (1-3), (1-4) o-nd 
(1.5): 

u = 27'tz'sin7 — 'j'w'^'^' cos 7 + £ + exp(e^|/sec7), (4.68) 

COS 2^ 

= £ tan 7 — 7' ■ro — 7'^^^' sin 7 + tan 7 exp(e^y sec7), (4.69) 

cos 7 

V = 27'e~^(cos7 tz^Q"' — sin 7) + ue'^ secy + 'y'{2uj'^' — tan'y)y (4.70) 

andp is given in (4-63), where (1)(4.48) is taken ande is given (4-46); (2) 7 is implicitly 
determined by (4-51), (4-52)-(4.54) are taken ande is given (4-46)- 

Next we consider the case (/g — z/x. Recall that we have already determined a in 
(4.27). Now (4.29) becomes 

x'^Tx — z'^Kz + Xz{Tz — Kx) — ZK -\- XT — 0. (4-71) 

equivalently, 

{xdx + zdz + 1) (r - -k) = 0. (4.72) 

Thus 

T = ^K + x~'^e{t,x/z) (4.73) 



for some two- variable function s. Denote 



^ = -. (4.74) 
z 
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Substituting (4.27) into (4.17), we get 

0:7^ + k{1 + W'^') — Tw'^'^' — x{l-ix + Tz) = 0. 

By (4.73), it changes to 
Thus 



w 



K = x'f'i + X w {e^' — H — . i^{t, w), 

vl + 



(4.75) 



(4.76) 



(4.77) 



where is another two- variable function and the last terms is so written for later conve- 
nience. Moreover, 



T = + X ^-uj{e'^' — £^) H — " -I- X ^e. 



1? 



For convenience of solving the problem, we denote 



So 



K = X'-^'i + X ^(f)-\- 



T = Z'-y'i + z '^ip + 



Observe that 



Kz — —X W (p^ — X ' 



2^3 / ' 



-2-^2, , -1 / ^^ti 



Thus 



and 



Tx~ >'^z — x'''^w^{;^ -|- 0)^ -|- x~^w\Jl -|- w"^ -Qi 



Moreover, 



(4.78) 

(4.79) 
(4.80) 

(4.81) 
(4.82) 

(4.83) 
(4.84) 



(4.85) 



and 



t(Tj; - Hz) ^ X W i))^ X 



+x-^ro[7i(V^ + 0)^ + + 7i\/l + w2 



(4.86) 
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Observe that 



( d 



(4.87) 



and 



x''^w^{wd^ + 2) 



(4.88) 



By (4.84), (4.87) and (4.88), (4.30) is equivalent to the following system of partial differ- 
ential equations: 



{wd^ - 1)(0(V' + = w^{wd^ + 3)(V'(V' + "0)^), 



(4.89) 



(^a^ - 1) (\/lT^0i?^) - ^"(t^a^ + 2) (x/r+^V^^?^) - , (V' + 0)^ = 0, (4.90) 



(4.91) 



+ 7i^^ = 0, (V' + 0)t^ = 0. 

According to (4.91), 

-d = a{t) + e-^Hi(ro), V + = /3i(^) + i(^), 
where Q;,/3i,ii and i are arbitrary one- variable functions. Write (4.89) as 

wd^\i{<^ — ^"^il^)] = o'{<j> + 'Uj'^ip)- 
Moreover, (4.90) is equivalent to 

w{l + w'^)d^[i\{(t) - w'^il^)] - L[{(f) + w^ip) - w^{ae^' + = 0. 
On the other hand, the first equation in (4.79) yields 

e = e^(^)£(i , w), = ^^e^(^)£^(i , w) , 
where £ is a function of t and w. According the second equation in (4.79) and (4.92), 

w''^(t) + w~^e = i3i + i^- e^[(l + w^Y^ + w'^^e] = p^ + t. (4.96) 

Thus 

(4.97) 



(4.92) 



(4.93) 



(4.94) 



(4.95) 



£ — 



W 



P + P: 



we 



V ^/(TT^ 
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dw 



wie 



7(1 + 



2^3 



dw 



for some function (3 oi t. Thus 



we 



^(1 + 



TJJ- 



2^3 



dw + 



^(1 + 



2^3 



Note 



7 V(l + ^2)3 7 7(1+^2)3 



P + P: 



I 



we 



dw + 



V v/(l + ^2)3 7 7(1+^2)3 



rote 



1 + ^2 ' 



(4.98) 



(4.99) 



t'(0 + ro^^) = .'((1 - ^2)0 + (A + Ot^') = ^ '^''"^^ 



1+^2 

— = (izi7 + / — — dw 



(4.100) 



Thus (4.93) forces us to take i! = 0. Replacing f3i by f3i + i, we have l = 0, and (4.93) 
naturaUy holds. Similarly, (4.94) yields Li = 0, and (4.94) naturally holds. 
According to (4.92), (4.95) and (4.98), 



,0 



— a, (j) = j3i — tl) — 



1+^2 



W^ 



P + Pl 



we 



dw 



(4.101) 



By (4.80), 



^2 



,0 



/we 
7?r+ 



1„9 



aro — Px e 



T = Z-f[ + piZ-^ ( . + 



1 + ^2 VI + ^2 y 7(1+^2)3 



/ 



we 



(izi7 + 



q; + Pz~'^e^ 



(4.102) 



(4.103) 



Modulo the transformation in (2.22) and (2.23), we take g — a ^. According to (4.18), 
(4.27) and (4.101), 



3,-3g-7i 



w^ 



P + Pl 



we 



'J ^y{l + w^f 



dw 



(4.104) 
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According (4.83), — t^. Hence (4.16) gives 

Recall ^ = f + Ky + ge^y and rj = ry + (e^K By (4.1), (4.27), (4.102) and (4.103), we 
have the following second theorem in this section. 

Theorem 4.2. Let w = x/z. Suppose that a, f3, (3i,'~fi are arbitrary functions oft and 
3" is any function ofw. We have the following solutions of the three-dimensional classical 
non-steady boundary layer equations (1-3), (1-4) o-nd (1-5): 



u 



+ ^^7=i£i^ + exp (xye^+-0 ' (4-106) 
V 1 + 



I e^ f we 

w — z^^+ Piz ^ — H — , / — , = dw 



' + -exp(xye^+^i) , (4.107) 



({x" ^z'')'^'y y e"^"^! \ ^ f we~^ , \ c,,,, 

, ^ - 4 V- /5 + A / , dw + i/xe^+Ti 



+ '^''7:^'°;-' - ( 2,; + , - i exp (.,e^-'.) (4.108) 

x(a;^ + 2;^) V ^/x^ + z^J x 



and p is given in (4-105) with k in (4-102) andr in (4-103)- 

5 3-D Distinct Exponential Approach 

In this section, we will find certain function-parameter exact solutions with two distinct 
exponential functions in y for the three-dimensional non-steady boundary layer equations 
(1.3)-(1.5). 

We use the settings in (4.1)-(4.12) and continue our case-by-case approach. 
Case 2. ^'^ = e^. 
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In this case, 

^y^yx - i^x + Vz)^yy + Vyiyz = nnx + - a({ag)z - a^gC^ - aQK^e'"^ + {ag{ag)^ 

-cr'^9gx)e^'^ + {{crgK,)x + {crg)zr + crg[{K,a^ + ra^ - (7{n^ + T^))y - (7/^]}e^ (5.1) 

by (4.11) and 

CyVyx ~ i^x + 'nz)Vyy + VyVyz 

= KT^ + TTz - (yg{(yC)x - (r^gxC + crgT^^e"^ + {aC{aC)^ - a^COe"^"^ 

+{-{(tCt)z - {(tQxK, + (tC,[{k.(7x + - (7(k^ + Tz))y - (Tfx]}e~'^ (5.2) 

by (4.12), where we treat w = ay. Moreover, (4.14) holds and 

Vyt = Tt + (aatCy - icrC)t)e~^, Vyyy = -cr^C^'^- (5-3) 
Thus (4.2) and (4.3) are imphed by the following system of partial differential equations: 

Kt + KK^ + TKz - aC{ag)z - a'^gCz + -Px = 0, (5.4) 

P 

Tt + KT^ + TT^ - agiaOx - cf'^QxC, + -Pz = 0, (5.5) 

P 

i^z = 0, = 0, (5.6) 

CTx = 0, az = 0, (5.7) 

at + K,ax + raz — a{K,x + Tz) = 0, (5.8) 

{ag)t + {agn):c + {'79)zr - cr^g{fx + i^cr) = 0, (5.9) 

{aOt + {(^Cr), + {aOxi^ + a\{U - = 0- (5.10) 

According to (5.7), 

a^^{t) (5.11) 
for some nonzero function 7 of t. Moreover, (5.8) can written as 



l^x + Tz^-. (5.12) 

7 



By (5.6) and (2.17)-(2.20), we take 



K = r = a 2;, (5.13) 

7 

where o; is a function of t. Modulo the transformation in (2.22) and (2.23), we take 
g — 7"-*^, and (5.9) becomes 

fx 1^1. (5.14) 

7 
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Substituting (5.11) and (5.14) into (5.10), we get 

2(7' - ^1^)C. + 7Ct + a'7<. + (V - a'7)<. = 0. (5.15) 

Thus 

C = ^0(e"x/7,e-"z). (5.16) 

7 

On the other hand, (5.4)-(5.6), (5.11) and (5.13) say that the compatibihty p^z — Pzx 
in (5.4) and (5.5) is imphed by 

(C).. = (C)..- (5.17) 
Let e = ±1. The above two equations imply 

7 = ee^°^, C = exp (2^ J e^'^dt - 4q;^ [9(e-"(x + z)) + i(e-°^(x - z))], (5.18) 

where and l are arbitrary one- variable functions. Moreover, (5.4) and (5.5) yield 

p= f{.e.p(2./e-*-2a.)p(e-"(...))-,e-"(.-.))l 

-{a" + {a')^){x^ + z^)} (5.19) 

modulo the transformation in (2.21). 

By (4.1), we have the following theorem. 

Theorem 5.1. Let l be one-variable functions and let a be function oft. Suppose 
e = ±1. We have the following solutions of the three-dimensional classical non-steady 
boundary layer equations (1.3), (I.4) o,nd (1.5): 

u^a'x + exp(ee^"y) , (5.20) 

» ^ - eexp (-ee-. + 2. / - 2a) P(e-«(. . .)) . - .))]. (5.21) 

V = e(i/e^" - a'e^") - 2a'y - exp (^-ee^"y + 2p j e^'^dt - bo^ 

X [^'{e-^ix + z)) - i\e-"{x - z))] (5.22) 

and p is given in (5.19). 

Case 3. = and $ = = C- 

As in the earlier cases, we take ag — 1. Then 

^y^yx " {^x + Vz)^yy + Vy^yz = ^^x + + '^x'^ ^ 

+{nx + [{ncr^ + - a{K^ + Tz))y - af^]}e'^. (5.23) 
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by (4.11) and 

^yVyx - {L + 'nz)'nyy + -qy-qy:, = + TT^ + agr^^e^. (5.24) 

by (4.12). Thus (4.2) and (4.3) are implied by the following system of partial differential 
equations: (4.16), 

(7a; = = Tr, Gt + Ta ^ - a{Kx + T^) = l^x = cr{fx + i^a). (5.25) 

We write 

a = e'^'''\ T = e{t,z) (5.26) 
for some functions q and £ of i and z. The second equation in (5.25) yields 

>^x = + — £z = '4'i't, z) is a function of t, z =^ k, = 4>{t, z) + ■^(t, z)x (5.27) 

for some function ^ of i and z. The compatibility p^z — Pzx in (4-16) is equivalent to 

dz{K,t + K'K'x + 'TK'z) — (5.28) 

Note 

KKx + TKz — (f)t + i^t^ + "00 + V'^^; + e(j)z + £'02^:. (5.29) 
For simplicity of solving (5.28) for 0, ^0 and (5.27) for c;, we assume 

0t + + £0z = 0, V't + V'' + sV'z = 0, £ = (5.30) 
for some functions a of i and t of ^. Thus 

V^= ,^ J ^ V '^ = 7X^^' <? = ln»2(« + 0(^ + ^(« + 0)-ln6' (5.31) 
for some one- variable functions 3^, Hence 

+ + r = 4, (5.32) 
^ '^2{(^ + i){t + %a + i)) ^ i[ 

-f(-'7f-(7)l (-) 

modulo the transformation in (2.21). Furthermore, the last equation in (5.25) yields 

Recall ^ = / + + ge!^^ and r] = ry. By (4.1), we obtain the following theorem: 
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Theorem 5.2. Let l be any function of z and let a be any function oft. Suppose that 
and ^2 a^^e arbitrary one-variable functions. We have the following solutions of the 
three-dimensional classical non-steady boundary layer equations (1.3), (1.4) o,nd (1.5): 

u = + + exp[y(.')-^^>2(« + i)(t + ^(a + l))], w = (5.36) 

t + ks{a + t) l' 

V^-<^Ja + i)(t + '^(a + t))-— rT^—7Z7 T^T 7 (5-37) 

and p is given in (5.34). 

Case 4- H = and $ = e'^^ for a function 7 of t. 
Again we assume ag — 1. In this case, 

+ [(ko-x + Tcr^ - a{K^ + T^))y - af^e"^ + a^ge^'^ . (5.38) 
by (4.11) and 

iy'Qyx - {ix + ri:,)riyy + r]yr]yz = HT^ + + crgT^ce"^ 

+7Va,C'e^^- + ia{g{aOx - 7<7^.C)e('+^)^. (5.39) 

by (4.12). Moreover, we have (4.14) and 

riyt = rt + {{laQt + ^aatCy)e'^^ , rjyyy = (7^)^(6^^. (5.40) 

Thus (4.2) and (4.3) are imphed by the following system of partial differential equations: 
(4.16), 

Hz = Cx = Cz = (^z = (^x = = 0, (5.41) 

CTt - (j(«x + Tz) = 0, a{fx + va) (5.42) 

{-laQt + 7aCr, = ^""a^if^ + u^a). (5.43) 
According to (5.41) and (5.42), 

cr = /3, fi; = a'a;, r = — — ^ (5.44) 

modulo the transformations in (2.17)-(2.20), where a and /3 are arbitrary function of t. 
Moreover, 
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Then (5.43) becomes 



(In7/?C)* + ^^-^ = l{a' + i/(7 - (5.46) 



(5.47) 



So 

/- - J_„a+/7(a'+i^(7-l)/3")rft 

Now 

p = -|[(a" + (a')')a^' + {(^" - 2«'/?' + ((« )' - «V)/^"''s'] (5-48) 
modulo the transformation in (2.21). Moreover, 

^^f + a'xy+^e^y, (5.49) 



^ _ yr- ^t^jy- _|_ ^a+/j7;/+.l 7(,a +!^(.7-l)/5"jat_ (5.50) 

By (4.1), we have the following theorem. 

Theorem 5.3. Let a, (3 and ^ be any functions oft. We have the following solutions 
of the three-dimensional classical non-steady boundary layer equations (1-3), (1-4) o-nd 
(1.5): ^ ^ 

u^a'x + v^v(3-'^-^y, (5.51) 



_ {13' - a'l3)z ^ ]^^„+^^j,+j^(„/+^(^_i)^2 
and p is given in (5.48). 

6 3-D Trigonometric and Hyperbolic Approach 

In this section, we will find certain function-parameter exact solutions with trigonometric 
and hyperbolic-type functions in y for the three-dimensional non-steady boundary layer 
equations (1.3)-(1.5). 

We start with the general approach given in (4.1)-(4.12) and continue the case-by-case 
process. Again we use the notion zu — a{t, x, z)y. For a e IR, we denote 

^0 = ^ , = ^ , (6.1) 

■j?i = sinro, i^i^cosro. (6.2) 



Case 5. H = 'd^ with r = 0, 1, and ( = = $. 



28 



Note 

-{-iyag{[Kaa: + rcr^ - (7(«;a; + r^)]?/ - (7/^}i?^ + cra^g^-dl (6.3) 
by (4.11) and 

^yVyx - {^x + Vz)Vyy + VyVyz = + TT^ + C^g'T^iy,. (6.4) 

by (4.12). Moreover, 

^yt = + (o-g')A - {-lyaatgy&r: iyyy = (-l)^^^'!?^, r)yt = Tt, 77j/2/J/ = (6-5) 

by (4.4) and (4.5). Thus (4.2) and (4.3) are implied by the following system of partial 
differential equations: 

Tx^CTx^ fx^ 0, (6.6) 

1 1 

nt + nnx + + (a5o,r + Si,r)cr^ggx + -Px = 0, n + tt^ + -pz = 0, (6.7) 

{ag)t + iagK% + {ag)zT - {-lYua^g = 0, (6.8) 

C^t + T(7z — cr{Kx + Tz) = 0. (6.9) 

We take / = by (6.6). Let a be a function of t and let t be a function of z. Denote 

w — a — w = —=. (6.10) 
In order to solve (6.9), we first assume 

r=^, ^=|^x. (6.11) 

Then 

a ^ (6.12) 
for some one-variable function $5. Multiplying (6.8) by 2(7g, we obtain 

[{agf]t + K[{ag)\ + T[{agf]z + 2{k, - {-Ifva^) {agf = 0. (6.13) 

So 

{agf = l,(z^, exp (^(-1)^2^^^ j ^) (6.14) 

for some two-variable function <;. By the compatibility Pxz — Pzx, dxdz{crgy — 0. Hence 
we take 



(^{-iyM^{zu)r J (6.13) 
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where </? is a one- variable function. According to (6.5) and (6.9), 



p / 2a'a" - {ay ^ {a'Y , f dz 



(0 



(6.14) 



modulo the transformation in (2.21). 

Recall ^ — ny + gH{w) and 77 = ry. By (4.1), we have the following theorem. 

Theorem 6.1. Let l he any function of z and let a he an arhitrary function of t. 
Suppose that '^s, </? are any functions of w — a — l and a e R. We have the following 
solutions of the three-dimensional classical non-steady houndary layer equations (1-3), 
(1.4) and (1.5): 

a at y a 



p is given in (6.14), and 
ii 

a if 

u — -TT—.X + 



2a' 2v^ 

exp 



(y 2a' 



dz 



X 



a'^y 



+ a exp 



'a'^y 



or 



u 



a 

2d 



-.X — 



exp — 2z/Q= 



02 



dz 

w 



COS 



a"^y 



(6.16) 



(6.17) 



Suppose 

K^a'(t)x, r^(3'(t)z, a = e°+^ (6.18) 

for some functions a and (3 of t. Equation (6.9) naturally holds. By the compatibility 
Pxz = Pzx in (6-7), we take 

{agf = e{t,z)-q{t,x) (6.19) 

for some two- variable functions e and <;. Now (6.13) is implied by the system of the 
following equations: 

£t+/3'^£, + 2(Q;'-(-l)Ve2(°+^))£ = 0, ^t + Q;'x^^ + 2(Q;'- (-l)Ve2(°+^))^ = 0. (6.20) 



Thus 



and 



e = '~s{ze-^) 



^ = t(xe-") 



,e.p(-2„ + (-ir2./e--.*) 
exp(-2„.(-l,-2./e--..,). 



(6.21) 
(6.22) 
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where t and $5 are arbitrary one- variable functions. Hence 



, ^ V.(.e-.)-4.e-.)exp - 2„ . (-1) V / e=<«-.*) 



(6.23) 



and 



-(a" + (a')V-(/9" + (/9')')^'] 



(6.24) 



modulo the transformation in (2.21). 

Recall ^ — Ky + gH{w) and rj — ry. By (4.1), we have the following theorem. 

Theorem 6.2. Let l, $5 be any one-variable functions and let a, j3 be any functions of 
t. We have the following solutions of the three-dimensional classical non-steady boundary 
layer equations (1.3), (1.4) CLnd (1.5): (1) w = P'z, 



u = a 



: + V3(..-)-^(^e-.).xp (-„-./ e^<-«*) o„s(e-.). (6.25) 



e-°'t'{xe-'') 



exp 



X sin(e"+^y) - (a' + /3')y 



and p is given in (6.24) ^^^^ ^ — 1; (^) = /^'^; 

X [exp(e"+^y) + aexp(-e"+^y)] 



(6.26) 



(6.27) 



e-"t'(xe-") 



.2(a+/3) 



(it 



, exp ( —3 — 2a-\- u / e 

4v/S^(ze-/^) - i(xe-") ^ ' 

X [exp(e"+^y) - a exp(-e"+^y)] - {a' + (3')y (6.28) 
for a e and p is given in (6.24) ^^^^ r = 0. 

Case ^. = $ = in (6.1) and (6.2), and a — e"*^*'' for some function a of t. 
Note = e"|/ in this case. Moreover, 

+e"[(^«:), + g,T + C«:.]^r + (a^o,. + 5i,r)e'''g{g, + G) + e'"(C5. - ^G)^' (6-29) 
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by (4.11) and 

= /^t + {a'g + gt)e'^^r + {-Ifa't^'^gy^r. iyyy = {-Ife^'^gK. (6.30) 

By symmetry, (4.2) and (4.3) are implied by the following system of partial differential 
equations: fx — 0, 

a' — Kx — Tz — 0, (6.31) 

Cgz - gCz = 0, gCx - QxC = o. (6.32) 

+ i^i^x + + (a5o,r + h,r)e^°'g{gx + G) + ^Px = 0, (6.33) 

n + /tra; + rr^ + (a5o,r + 5i,r)e^"C(fl'a: + Q + = 0, (6.34) 

o^'g + gt + {gi^)x + g^r + (i^z - (-i)Ve'"^ = O, (6.35) 

«'C + Ct + {Cr)z + + - (-l)Ve2"C = 0. (6.36) 
We take / = by (4.1). Moreover, (6.31) yields 

a'x - T = (6.37) 

for some function <; of t, x, z. By (6.32), 

C, = (7 tan 7 (6.38) 

for some function 7 of t. According to (4.28) and (4.29), (6.35) and (6.36) give 

7' sec^ J + Tx — Kz tan^ 7 + (r^ — Kx) tan 7 = 0, (6.39) 

equivalently, 

7' sec^ 7 + <5"a;a; + <b"00 tau^ 7 + 2^^^ tan 7 — a' tan 7 = 0. (6.40) 
Now we use the notations in (4.31). The above equation is equivalent to 

d^{<i) — a' sin 7 cos 7 — 7', (6-41) 

that is, 

— (f){t, -uj) + ilj{t, -uj)w -\ — (a' sin 7 COS7 — 7')g7^ (6-42) 
for some functions and ip oi t and w. So 

K — a'x — {(f)^ + wip.^) cos 7 — sin 7 — {a' sin 7 cos 7 — 7')^ sin 7, (6.43) 

r = — (0^x7 + tiJ'^^) sin7 + C0S7 + (a' sin7 C0S7 — 7')tz7cos7. (6.44) 
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Note 

Tx - i^z^ (t>^w + '^i'^vj + Oi sin 7 cos 7 - 7'. (6.45) 
For simplicity of solving the problem, we assume 

^ = -^(Q;'sin7 COS7 - 7')ti7^, V'txT^zr = 0, (6.46) 

that is, Tx — Kz — 0. The compatibility p^z — Pzx in (6.33) and (6.36) is equivalent to 
d^d^ig^) = 0. So 

^e{t,uj)+q{t,w) (6.47) 

for some function e of t and w, and some function ^ of t and ui. On the other hand, (6.35) 
can be written as 

gt + i^gx + rgz + {a' + Kx + Kz tan 7 - (-1)^6^")^ = 0, (6.48) 

equivalently, 

gt + a'xdxig) - {(j)^ + '^i^^)d^{g) + (^ + (a' sin 7 C0S7 - i)w)d^{g) 

+(q;'(2 - sin^ 7) - + 7' tan 7 - (-l)Ve2")c/ = 0. (6.49) 

By (4.31) and (4.32), we obtain 

gt + [w{a' cos^ 7 - tl)^) + 7'ro]9^(g') + {a'w s\v? 7 - 7'ro + ilj)d^{g) 

+ {a{2 - sin^ 7) - + 7' tan 7 - (-l)Ve2")c/ = 0. (6.50) 

Multiplying 2g to the above equation, we have 

{g^)t + [w{a' cos^ 7 - i;^) + y-cu] (51^)^ + {a'w sin^ 7 - -f'w + ^) (g^)^ 
+2{a\2 - sin^ 7) - -0^ + 7' tan 7 - (-1)^6^")/ = 0. (6.51) 

Substituting (6.47) into the above equation, we get 

£t + + [^(a' cos^ 7 - V'tzr) + 27'^]^^ + {a'w sin^ 7 - 2'y'w + V')£tS7 

+2(a'(2 - sin^ 7) - + 7' tan 7 - (-l)Ve2")(£ + ^) = 0. (6.52) 

In order to solve (6.52), we assume 

e^^2pw + pi, q^^2l3w + l32 (6.53) 

for some functions of So £ is a quadratic polynomial in w and is a quadratic 

polynomial in tu, whose coefficients are functions of t. Now (6.52) is implied by the 
following system of partial differential equations: 
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et + 2/^27'^ + («'^ sin^ 7 + V') C^P^ + Pi) 

+2{a'{2 - sin^ 7) - -0^ + 7' tan 7 - (-1) Ve^")^ = 0, 

- 2/3i7'^ + [^(a'cos' 7 - Vt^)](2/3^ + P2) 
+2(a (2 - sin^ 7) - t/,^ + 7' tan 7 - (-l)Ve2")<? = 0. 

By the coefficients of quadratic terms, we have 

/3' + 2(2a' + 7'tan7 - (-IjVe^'^)/? = 0, 

/3' + 2(a'(2 + cos27) -2V'^ + 7'tan7- (-l)Ve2")/? = 0, 

which implies 



For simplicity, we take 
According (6.56), 



, a 

W€j^ -i^ cos 27. 



■0 = —07 cos 27. 



exp 



-lY2v J ( 



cos^ 7 



for bi e R. 

The coefficients of w in (6.54) and of w in (6.55) give 

+ 2/327 + (7a;72 + 27' tan 7 - (-l)''2i/e'")/3i = 0, 
- 2(3ii + (7a72 + 27' tan 7 + 2ue^'')i32 = 0. 



So 



/3i = (62 cos 27 + 63 sin 27) 



P2 = {-h sin 27 + 63 cos 27) 



exp ({-iy2iy J e^^'dt - 7a/2 
exp (^i-lY2iy J e^'^dt - la/2^ 



cos 7, 



cos^ 7 



(6.54) 

(6.55) 

(6.56) 
(6.57) 

(6.58) 

(6.59) 

(6.60) 



(6.61) 
(6.62) 

(6.63) 
(6.64) 



for 62, ^3 £ R- According to (6.47), we take the zero constant term of ^ and the constant 
term of e to be 



exp [{-iy2v j e^'^dt - 3a 



cos^ 7 



(6.65) 



by (6.54) with 64 e M. Hence we take 



g = ^J bie-"{w'^ + w"^) + e-°'/'^[{b2W + bsw) cos 27 + {hw - 1)2^) sin 27] + 64 



X 



exp n-l)V j e^'^dt -3a/2 



cos 7. 



(6.66) 
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Furthermore, 



K — — ^'{z cos 27 — X sin 27), 

Ot^ z 

T — — 7'(x COS 27 + 2^ sin 27) 



(6.67) 
(6.68) 



by (4.31), (6.43), (6.44), (6.46) and (6.58). 
Set 



Q^' / 9 9x 7' sin 27 ,99, / ^ 

K — (a;^ + 2;^) + - — - — -{x^ - z^) - xz^' cos2^. 



Then 

Moreover, 

e'"^(^. + G) 



(6.69) 
(6.70) 



cos 7 



[26ie-"ro + e-''/^(63 cos 27 - 62 sin 27)] 



exp 



(6.71) 



By (6.33), (6.34) and (6.69)-(6.71), 



r / /' o , / / o 0^^2-0^, 7" sin 27 + 2(7')^ cos 7 2 2^ 

p = p|a;z(7 cos 27 + a 7 cos 27 — 2(7 ) sm27) H [z —x) 



I' 



exp ( {-iy2u I e'^'^dt - a 



X [616-"^^ + e-"/2(63 cos 27 - 62 sin 27)^]} 



(6.72) 



modulo the transformation in (2.21). 

Recall C, = Ky + gH{w) and t] = ry + C,^ . Moreover, + w'^ = x"^ + z^ according to 
(4.31). By (4.1), we have the following theorem. 

Theorem 6.3. Let a, 7 he any functions oft and let a, 61, 62, ^3, ^'4 £ I^- For r = 0, 1, 

we define anddr in (6.1) and (6.2) withw = e"-y, andw andw in (4-31). We have the 
following solutions of the three-dimensional classical non-steady boundary layer equations 
(1.3), (1.4) and (1.5): 

0/ X I 

u ^ -1^ + \J &ie-"(a;2 + z"^) + e-"/2[(62^ + 63^) cos 27 + (63^ - h2w) sin 27] + 64 



X 



exp 



e^'^dt -a/2 



COS7 — 7'(2;cos27 — xsin 27), 



(6.73) 



w 



'z I : : 

— + W 6ie-"(x2 + z"^) + e-"/2 [(621^7 + 63^) cos 27 + (63^7 - 62^7) sin 27] + 64 



q; 2; 

X 



exp (-l)V / e"^dt-a/2 



sin 7 — 7' (a; cos 27 + 2; sin 27) , 



(6.74) 
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2hic-"^ + c-"/-(h:^ cos 2- - ho sin 2- ) 



2-s/bie-°'{x'^ + z^) + e-°/2[(62^ + 63^) cos 27 + (63^ - 62^) sin 27] + 64 



X 



X 



1?^ - a'y (6.75) 



exp y" e^'^dt - 3q;/2^ 

anc?p is given in (6.72). 

Case 7. H = dr and $ = -i?^ in (6.1) and (6.2) with r = 0, 1. 
In this case, 

CyCyx i^x ~l~ Vz)^yy ~l~ Vy^yz 

= + TK^ + (a(5o,r + 5i^r)(y'^gg^ + [(o-gi/t)^ + ((T5()^r]'(9^ + aa^g'^^l + (-l)''cr{Cfi:2 

+5'[('«o-a; + ra^ - a{K^ + r^))y - af^]}'&r + (-l)''K(f^5')2 - cr'^9Cz]'&r'^r (6.76) 

by (4.11), and 

^yVyx i^x ~l~ Vz^Vyy ~l~ VyVyz 
= KTa; + rr^ + (-1)''{[((tCi")2 + (o'C)2:/«]^r - {aSo,r + 5l,r)0-\Cj + ^-^-^C^^r + (^{O'^x 

+ {-iyC[{'^a, + ra, - a{K, + T,))y - aQ + (-l)'^[^(aC). - a^.C]^r}^r (6.77) 
by (4.12). Moreover, 

Cj/t = «t + (cr5')t1?r + {-lyaatgy&r: ^yyy = (-l)V^5ri?r, (6.78) 

?7,t = n + {-lYiaOt^r + {-lyaatCy^r, Vyyy = (^'C^r (6.79) 

by (4.4) and (4.5). Thus (4.2) and (4.3) are imphed by the following system of partial 
differential equations: 

CTx^cTz^ 0, - agfa; = 0, gr^ - (-l)'"(7C/x = 0, (6.80) 

Kt + KKx + TKz + (a(5o,r + Si^r)cr'^g9x + ^Px = 0, (6.81) 

Tt + KTx + TT, - {-iy{aSo,r + Si,r)a^CCz + -Pz = 0, (6.82) 

{ag)t + {<jgK)^ + {ag),T - {-lYva^g = 0, (6.83) 

{aQt + Ht), + (aC)./^ - (-l)Va\ = 0, (6.84) 

(7t - cr{Kx + T^) = 0, (6.85) 

C,9z - 9Cz = 0, gCx - 9xC = 0. (6.86) 
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By (6.80) and (6.86), we have 

cr = e", ( = gtanj, r^; cot7 — (— l)''fi;2tan7 = 0. 
Moreover, (6.83) and (6.84) become 

gt + i^9x + rg, + {a' + - (-1)^6^")^ = 0, 
Vg^sec^ -f + {gt + HQx + rg^ + {a' + - (-l)Ve^")5() tan7 = 0. 
Furthermore, [(6.89) — tan 7 x (6.88)]/^' yields 

7' sec^ 7 + (tz — f^x) tan 7 = =^ ~ = 27' esc 27. 
On the other hand, (6.85) says 



Thus 



We take 



Kx + r^^a. 



= — + 7 CSC 27, = — - 7 CSC 27. 



K 



a' 



+ 7' CSC 27 ) X, = 7' CSC 27 I 2; 



a 



(6.87) 

(6.88) 
(6.89) 

(6.90) 

(6.91) 
(6.92) 
(6.93) 



for simphcity. So the third equation in (6.87) naturally holds. Moreover, the compatibility 
Pxz = Pzx in (6.81) and (6.82) is imphed by {g'^)xz = 0. So 



g"" ^(t>{t,x)+i^{t,z) 
for some two- variable functions (j) and ■0. By (6.93), (6.88) becomes 



(6.94) 



a 



)i + TT + V CSC 27 x{g )a; + 17 - 7 CSC 27 z{g ) 



a 



Thus 



+ (3a' + 27' CSC 27 - {-lY2ue"^)g' = 0. 



g2 ^ e-^''+^-^^^^''^''"'^\i{xe-''/^y/^)+^{ze-'''^y/i^))coi-i 



(6.95) 



(6.96) 



for some one-variable functions t, According (6.81), (6.82), (6.87), (6.93) and (6.96), 
we have: 

p = ^{(a5o,^ + 5i,^)[(-l)^(tan7)$5(ze-"/^Vta^n') - (cot7)i(xe-"/^Vcot7)] 



^g-a+(-l)'"2i//e2«dt _ 



/ \ 2 // 

^ + 7' CSC 27 ) + ^ + CSC 27 (7" - 2(7')^ cot 27) 



X 



— - 7' CSC 27 ) + ^ - CSC 27 (7" - 2(7')^ cot 27) 



(6.97) 
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modulo the transformation in (2.21). Furthermore, we take / = by (6.80). 

Note ^ — Ky-\- g'&r and rj — ry + By (4.1), we have the following theorem. 

Theorem 6.4. Let l, ^ be any one-variable functions and let a, 7 be any functions of 
t. For r — 0,1, we define and -dr (6-1) and (6.2) with w — and a e R. We 
have the following solutions of the three-dimensional classical non-steady boundary layer 
equations (1-3), (1-4) o^nd (1.5): 

u = (^^+ycsc27^x + e-°/2+(-^)^'^/^'"'^*i 

x^(t(a;e-°/2Vcot7) + a(^e-"/2Vtair7))cot7, (6.98) 

w = (^^-ycsc27^z + (-l)^e-"/2+(-ir-/e^'^<^*^^ 

x^(t(xe-"/2V^^) + ^(^e-"/2Vtan7))tan7, (6.99) 

„-2a+(-l)''!^/e2«dt 

V — — q; y — 

2^(i(a;e-"/2^/5ot7) + Q{ze-'^/^^/ta^l^)) 

x[(cot7)i?^i'(a;e-"/2v/cot7) + {tan-f)^rQ'{ze-''/^^/ta^)] (6.100) 
and p is given in (6.97). 

7 3-D Rational Approach 

In this section, we will find certain function-parameter exact solutions of rational in y 
for the three-dimensional classical non-steady boundary layer equations (1.3), (1.4) and 
(1.5). 

First we consider another three cases of solving the three-dimensional classical non- 
steady boundary layer equations (1.3), (1.4) and (1.5) based on the settings in (4.1)-(4.12). 



Case 8. H = ^ = y-\ 

In this case a — 1 and uj — y. So (4.11) becomes 

-2g[{K^ + T,)y + f^]y-^ + (Cy, - gg^ - 2gQy-\ (7.1) 
Moreover, 

Cj/t ^n- gty~^, Cyyy = -Qgy'^- (7.2) 
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By symmetry, (4.2) and (4.3) are implied by the following system of partial differential 
equations: fx — 0, 

Kt + KKx + TKz + -Px = 0, Tj + KT^ + TT^ + = 0, (7.3) 

P P 

C9z - 99x - 2^C. + 6i^^ = 0, gCx - CCz - KQx + = 0, (7.4) 

9t + C«2 + f^9x + 9zT + 5'(3Ka; + St^) = 0, (7.5) 

Ct + 9rx + «;Cx + rC. + C(3r, + 2«;^) = 0. (7.6) 

We take / = and assume C, — g'tan7 for some function 7 of t. Use the notations in 
(4.31). Then both equations in (7.4) are equivalent to the equation: 

9x + 9z tan 7 = 6z/ =^ g — Gux + (j){t, w) (7.7) 

for some two variable function 0. Moreover, (7.6) — tan 7 x (7.5) is implied by 

7' sec^ 7 + Tj; — ^2 tan^ 7 + (t^ — Kx) tan 7 = 0, (7.8) 

equivalently, 

7'sec7 + 9^(t — «;tan7) = (7.9) 

by (4.32). So 

T = fi;tan7 + ■i/'(i,ti7) — 7'^sec7 C^-IO) 
for another two- variable function ■0. To simplify the problem, we assume 

«z - = ~ 9^(k) sec 7 = -tl)^ sin 7 - 7' (7-11) 

by (4.31) and (4.32). We take 

K = aro — -■0 sin 27 — 7'ro cos 7 (7-12) 

for some function a oit. 
Note that (7.5) becomes 

0t — ^'w^^ + 6z/Q;tz7 — 3i/'0 sin 27 — ^v^'w cos 7 + (-0 — ^'w sec 7)^^ cos 7 
-|-(6i/x + 0)(3Q;sec7 + 2'0^cos7 — 27'tan7) =0, (7-13) 

equivalently, 

<^t — Si^V' sin 27 + ((/) — 6z^'iX7 sin 7) (3q; sec 7 + 2'^^ cos 7 — 27' tan 7) — ^v^'w cos 7 
-|-0^'0cos7 + 2ot(12i/q; — 7Vro + 6i/'0^cos^7 — 61/7' sin 7) = 0. (7-14) 
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First we have: 

12i/Q; — 7Vra7 + 6i^'0ro cos^ 7 — 61/7' sin 7 = 0. 

So we take 

th — —7V sec^ 7 + 7'tz7 tan 7 sec 7 — 2q;g7 sec^ 7. 

Now (7.14) becomes 

(ht + (cos 7 — 3i/ sin 27) ( — 7'^ sec^ 7 + 7'tz7 tan 7 sec 7 — 2Q;n7 sec^ ' 

—Giyj'w cosj H 7'0^(0 — 6i/-±Jsin7) sec 7 = 0, 

equivalently. 

0t + sec 7 — 7'^ tan 7 — (7' tan 7 + 20; sec ^)w(f)^ 

+6i/(2Q;tan7 — 7'sec7)ro = 0. 

For simplicity, we take 

for a function /3 7^ 6 sin 7 of Then (7.18) gives 

+ 2i/(/?' cos 7 - 2/5y sin 7 - 67') 
4i/(/? — 6i/sin7) 

In summary, we have 

g — ^vx + C = (61/x + /3ro) tan 7, 

K = aw — ( — 7'/? tan 7 + 7' sin 7 tan 7 — 2q; tan 7 1 — 7'-^' cos 7 
\6u J 

/3^i + 2u{/3' cos 7 - 2l3i sin 7 - 6y ) ^ , , ^ . ^ 

= 7-:; ^ ro + 2rotan7 

4i/(/3-6z/sin7) ^ 

— f — /37'tan7 + 7'sec7 I ro, 
\6u J 

+ 2u{(3' cos 7 - 2/3y sin 7 - 67') , 

T = —7- — - — : r (ro tan 7 — 2ro) 

4i/(/3 - 6z/sin7) ^ ' ' 

H (3'-^''nj — ^'w sec 7. 

6i/ 



Set 



(j-'i + 2u{(3' cos 7 - 2(3i sin 7 - 6y ) , . , 

K = -t-t: : r [W —2vd) 

8i/cos7 (p — 6i/sm7) 

1 7' . 
H P"-)' sec 7 H tan 7('^^ — ro^) — ^'ww. 
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Then 

K — kx, T — kz- (7-25) 
By the fact Hz = t^, (7.11) and (7.24), we have 

p^-pkt-^{^ + T^) (7.26) 

modulo the transformation in (2.21). 

In this case, ^ = ny + gy~^ and i] = ry + Cy~^- By (4.1), we have the following 
theorem. 

Theorem 7.1. Let /3,7 be any functions oft. In terms of the notations in (4-31), we 
have the following solutions of the three-dimensional classical non-steady boundary layer 
equations (1.3), (1.4) o,nd (1.5): 

(3^y + 2u{(3' cos 7 - 2/?7^ sin 7 - Gj) ^ , , ^ . ^ 

u — : — 7- r (ro + 2rotan7) 

4z/(/?-6i/sin7) ^ 

- (^-^Pl' tan 7 + 7' sec 7^ - {6iyx + /3w)y~'^, (7.27) 

(3^i + 2t/(/3' cos 7 - 2/37' sin 7 - 67') _ ^ ^ . , 

4z/(/3-6i/sin7) ^ ' ' 

-{-^Pj'w — j'wsec^ — {61/ X -\- /9t£r)t/~^tan7, (7.28) 

/3V + 2i^(/5'cos7-2/37'sin7-67') ^ . 

V = — J . ' ' -y - sec 7 - 6uy-' 7.29 

4i/(p — bu sm 7) cos 7 bu 

and p is given in (7.26) via (7.22)-(7.24). 
Case 9. // = $ = 0. 

In this case, / can be any function of t, x, z, and (4.2) and (4.3) are equivalent to the 
equations in (7.3), whose compatibility p^z — Pzx is equivalent to (4.30). The simplest 
solutions are 

K — 9x{t,x,z), T — 9z{t,x,z) (7.30) 
for some three- variable function 9. Since Kz = Tx, 

P^-p9t-P-{9l + 9l) (7.31) 

modulo the transformation in (2.21). Besides, the k and r in Cases 1,3,4 and 8 work for 
this case. 
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Next we assume 



K — e'^°'z + hx{t,x, z), T — hz{t,x,z), 



(7.32) 



where a is a nonzero function of t and /i is a function of t, x, z. Now (4.30) becomes 



lOL + hxx + /izz = 0. 

Hence 

h — Tj^ — olx^ 

for some time-dependant harmonic function T(t,x, ^), that is, 



In this subcase, 



>^ — ^xx — '^ot'x + e^"^, T = Ta 
i^t ~ ^xxt — 2q;''x + 2a'e^"z, 

Kx '^xxx 2q; , xxz ~l~ 6 



2a 



(7.33) 
(7.34) 

(7.35) 

(7.36) 
(7.37) 
(7.38) 



Thus 



and 



So 



Kt + + TK^ = Txxt — 'ioi' X + 2a'e'^°'z 

+ {Txx - 2a'a; + e'"^)(T,,, - 2a') + T,,(T,,, + e^") 

T.t. + ^(T^, + T^J, - 2a' (xT,,), + A{a'fx + e2"(zT,,, + T,,) - 2a" x 

Tf -\- KTx ~\~ TTz — ^ xzt ~l~ xx 2q! X -\- C Z^^xxz ~I~ ^ xz'^ xzz- 



'^xt + 2(TL + TL) - 2q;'xT,, + {2{a'f - ol')x^ + e^'^(zT,, + T,) 



(7.39) 



(7.40) 



(7.41) 



modulo the transformation in (2.21) 
Denote 



2 I 2 

W — X -\- Z . 



K = a'x + z4>{t, zu), T = a'z — x(j){t, w) 
for some functions a of i and (f) oi t and zu. Note 

l^z-Tx^ 2{w(p)^. 



(7.42) 
(7.43) 

(7.44) 
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Then (4.30) becomes 

{zu(P)^t + 2a'[(ti7</.)^ + zu{w(f>)^^ = 0. (7.45) 

Thus 

zu 

for a function (5 olt and a one- variable function 9. So 

, z{^{e-'"^w) + (5) , x($>(e-2"ro) + /3) ,^ 

K = olx^ -^-^ ^ — T = a'z ^-^ — —. (7.47) 

zu zu 



Moreover, 

} x -r 

zu zu 



, ,1 / (^J(e-2"ti7) + /?)2 . 
Kt + + TK, = (a" + (a') )a; + ^ - — (^-^S) 



/ // / A2N (»(e-2"OT) + /?)2 . ^, 

Tt + + TT, = (a" + (q;')')-^ - ^ - -2^^ ■ (7.49) 



,'N2 

Hence 



2 t 7 [e-2"(a;2 + ^2)]2 

arctan - (7.50) 

X 

modulo the transformation in (2.21). 

The final subcase we are concerned with is to assume 

T^e{t,z) (7.51) 

for some two-variable function z. In this subcase, (4.30) is equivalent to 

52('«t + /^/^a: + TK2) = 0. (7-52) 

We look for a solution of the form 

= 4>{t,z)^i}{t,z)x (7.53) 
for some two- variable function and V'- Note 

Kt KKa; + = 0t + "0*2^ "00 + V'^a; £^02 Szi^zX- (7-54) 

So (7.48) is equivalent to the following system of partial differential equations: 

dMt + + £0.) = 0, 9,(V't + V'' + e^z) = 0. (7.55) 
To solve the above system, we assume 

£ = ^ - ^ (7.56) 
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for some functions a of and q oit and z. We have the following solution of the above 
second equation: 



1 

for another one-variable function 3^. By the first equation, 



^ - TTi:^ (7-57) 



for another one- variable function and a real constant h. In this subcase, 
By (7.3) and (7.54), 



dz 



modulo the transformation in (2.21). 

Recall i = f + ny and rj = ry. By (4.1), we have the following theorem. 

Theorem 7.2. Let /i, 9 he arbitrary function oft, x, z and let ^, be any one-variable 
functions. Suppose that a, (5 are any functions oft, q is any function oft and z, h is a real 
constant and T(t, x, z) is a time- dependent harmonic function in x, z. Then we have the 
following solution of the three-dimensional classical non-steady boundary layer equations 
(1.3), (1.4) and (1.5): (1) 

U^9^, W^9^, V ^ -{9^a: + Ozz)y + 1^ (7.61) 
and p is given in (7.31); (2) 

u = ^xx — 2Q;'a; -|- e'^°'z, w = T^z, v = iJ,-\- 2a'y (7.62) 
and p is given in (7.41); (3) 



,^ , z{'^{e-^ ^{x^ + z^)) + 13) 
x'^ -|- ^' 

,^ x(^(e- ^°(x^ + z^)) + /3) 
-|- z 

and p is given (7.50); (4) 



= v^^-2a'y, (7.63) 



w^a'z- "^"^^ (7.64) 



" = TTW) ■ "^-^^^W ? 
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and]) is given in (7.60). 
Case 10. H = ^ = y\ 
In this case, a — 1. Now (4.11) becomes 

^y^yx ~ i^x + Vz)^yy + Vy^yz — l^t ''f ^^^^x + "^^^z ~ "^Qfx 

+ 2[(5'fi;)a; + QzT + (K-z " gi^x + Tz)]l/ + 2(5'5'a: + '^CQz " Kz)?/^- (7-67) 

Moreover, 

Cyt = ^1 + "^gtV, iyyy = 0. (7.68) 

Modulo the transformation (2.22) and (2.23), we assume r = 0. By symmetry, (4.2) and 
(4.3) are equivalent to the following system of partial differential equations, 

1 1 
Kt + KK^ - 2gf^ + -p^ = 0, -2Cfx + -Pz = 0, (7.69) 
P P 

9t + i^9x + C«2 = 0, Ct + Cx« - CK'x = 0, (7.70) 
99x + X9z - 9Cz = 0, CC. + 2Cx^ - 9xC = 0. (7.71) 
After some computations, we take 

9^a{t), C-Pit) (7.72) 
for functions a, (3 of t. So (7.71) naturally holds. By (7.70), we take 

. = ^^. (7.73) 
The compatibility p^z = Pzx in (7.69) gives 

2{pd,-ad,)if,)^j. (7.74) 



Thus 



„ a'' {(3x — az) , , ^ ^ „ ,„ __s 

2/j(a2 + /32-) ^<Pvj\t,^), w^ax + l3z (7.75) 



for some two-variable function (b. Thus 



^ , , ^ , off' n pa''{ax^ + 2Bxz — az^) ,„ 
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modulo the transformation in (2.21). 

Recall i — f-\-Ky-\- gy^ and 77 = (y"^. By (4.1), we have the following theorem. 



Theorem 7.3. Let a, (3 be arbitrary function oft and let (f){t,'uj) be an arbitrary two- 
variable function. Then we have the following solution of the three-dimensional classical 
non-steady boundary layer equations (1-3), (1-4) o.^d (1-5): 

3'x — a'z a''(az — 6x) , , ^, d'y 

„ = 2„,+ ^L_, ^^2fy. „ = _i-j_£J-^,(,,^)-^ (7.77) 

and p is given in (7.76), where w — ax -\- (5z. 

By the transformation in (2.22) and (2.23), we next consider the solutions of (4.2) and 
(4.3) in the following form: 

^^fy'^ + gy + h, tj = ny"^ + Ty^ + ijy, (7.78) 

where /, g, h, k, r, fj, are functions of t, x, z with / 7^ 0. Note 

= + g, = 6/y, ^yyy = 6/, ^yt = 3/^^^ + , (7.79) 

^yx = + g^, ^y^ = Sf^y'^ + g^, = f^y^ + g^y + h^, (7.80) 

r]y = 3Ky'^ + 2ry + rjyyy ^ 6k, rjyt = SKty'^ + 2Tty + i^t (7-81) 

Vyx = 3Ka;?/^ + 2Ta;y + yu^, r]y^ = SK^y"^ + 2T^y + jj,^, r]^ = K^y"^ + r^y'^ + fi^y. (7.82) 
So we have 

CyCyx ~ i^x + Vz)Cyy + VyCyz 

= (3/?/^ + g) i^fxy'^ + gx) - Gfy{{fx + K'z)y^ + r^y^ + {g^ + Hz)y + /ix) 

+ (3/€i/2 + 2rt/ + ^)(3/,|/2 + (7,) 
= 3(//a; - 2/k^ + 3f^K)y^ + 6(r/^ - fr^)y^ + 3{f^g - fg^ - 2///^ + Kg^ + f^ii)y^ 

+ {2Tg^-6fh^)y + lig^ + gg^, (7.83) 

^yVyx i^x ~l~ Vz^Vyy ~^ VyVyz 
= (3/y^ + 5-) (3/^0;?/^ + 2T^y + ii:o)-{ {fx + Hz)y^ + T^t/^ + {gx + IJ'z)y + K) 

X (6«y + 2t) + (3«;y' + 2Ty + //) (3«;,y2 + 2T,y + //,) 
= 2(3/Ta; + 2k^t - fxr)y^ + 2{gTx + r^fi - rg^ - 3Khx)y + gfix + ^-^-z - 2hxT 

+3(3/ - 2fxH + nKz)y'^ + [3(//^x + f^xg + ^^^/^ - '^t^gx - f^l^z) + '^rTz]y'^. (7.84) 
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Thus the equations (4.2) and (4.3) are imphed by the following system of partial differ- 
ential equations: 

//^ - 2fK^ + 3f^K = 0, - 2f^K + KK^ = 0, (7.85) 

rfz - fr, = 0, + 2k,t - Ur = 0, (7.86) 

ft + fx9 - f9x - 2fiiz + i^9z + fzl^ = 0, (7.87) 

3(«;t + ffJ'x + 1^x9 + i^zlJ' - 2n9x - Kji^) + 2tTz = 0, (7.88) 

- 3fhx = 0, Tt + 9Tx + T^n - - SkH^ = 0, (7.89) 

1 1 

9t + IJ'9z + 99x + -Px = Qi^f, l^t + 9l^x + A^Ai^ - 2/ia;r + -p^ = Guk. (7.90) 
P P 

Let 7 be a function of t. Again we use the notations w — 2; cos 7 — xsin7 and 
w — 2; sin 7 + X cos 7. We take the following solutions of (7.85): 

f — (f){t,w) cos ^, K — (f){t,w) sin ^ (7-91) 

for a two-variable function (p. Since fx + Kz — 0, (7.86) is implied by 

T^a(t){t,w) (7.92) 

for a function a of t. Now (7.87) and (7.88) become 

-7'^ sin 7 [0t 0^(// cos 7 - sin 7) - - 20//^] cos 7 ^g^;, sin 7 = 0, (7.93) 

3[(7' + A*x)0cos7 (0t 0^(/xcos7 - sin 7) - 20g^ - 0//^) sin 7] 2tTz = 0. (7.94) 
Moreover cos 7 x (7.93) — 3 sin 7 x (7.94) yields 

3(7' -I- jix cos^ 1 — 9z sin^ 7 -I- {jiz — 9x) sin 7 cos 7] -|- 2q;^0^ cos^ 7 = 0, (7.95) 

equivalently, 

2 

9^(//cos7 - g sin 7) = --q;^0^cos^7 - 7'. (7.96) 
According to (7.91) and (7.92), the first equation in (7.89) gives 

hx^ :^9zsec-f. (7.97) 

Substituting it into the second equation in (7.89), we get 

a'(j) + a[(j)t + 0^(Aicos7 - g sin 7) - (j)d^{g) sec 7] = 0. (7.98) 
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Furthermore, (7.93) can be written as 

— 7'0sin7+ [0t + 0^(/xcos7 — g'sin7)] cos^ — (f)d^{g) — 20^2 (/x cos 7 — gi sin 7) = 0. (7.99) 
The above two equations imphes 

92(cos7 II — sinj g) = — - ( — cos 7 + 7' sin 7 j . (7.100) 



So 



z I a' 



//cos 7 — g'sin7 = — - ( — COS7 + 7' sin 7 ) + e{t^ x) (7.101) 
2 V a 



for some two- variable function e. Substituting it into (7.96), we obtain 
sin7 f a' , .A , 2 ^ ^ 



cos 7 + 7' sin 7 + cos 7 = — q;0^cos 7 — 7', (7.102) 



2 \a ' ' 'y - ' 3 
which shows that both and (j)^ are purely a function of t. Thus we take 

(t)^(3w (7.103) 

for a function P oi t. Hence we have 

//COS 7 — gism7 = — -a pcos 7 + 7 I + I — — sm27 H tan 7 ] w (7.104) 

\3 /\3 2 2a J 

modulo the transformations in (2.17)-(2.20). 

In order to solve (7.90) and (7.98), we assume a — 2 and 

^a'/3cos'7 + 2y = ^ /3 = - ^ . (7.105) 
3 4 COS"' 7 

Then 

/fC0S7 — gsm.^ = ^'w =^ /* = 5ftan7 + 7'tz7sec7 (7.106) 
by (7.104). Moreover, (7.98) becomes 

Thus ^ 

g= cos-^ + 2i sin^w + %l){t,w) (7.108) 

for some two-variable function ■0. The compatibility p^z — Pzx in (7.90) becomes 

ifJ'x - 9z)t + {gifJ'x - gz))x + {iJ'ilJ'x - gz))z + 2-i'{wgz)x sec^ 7 = -6i//3 (7.109) 
by (7.97) and (7.105). 
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Observe that (7.106) and (7.108) yield 

fifz = —sin 7 COS 7 + 27' sin 7 + V'^ COS 7, ji^ - Qz ^ ^' - 'ipu^sec-f. (7.110) 
7 

So 

/ Sill ^ 

{l^x - gz)t = 7" - 7' — 5— + i^i^w^ sec 7 - V'^t sec 7, (7-111) 
cos^7 

{gilJ'x- 9z))x + {lJ'{fJ'x- 9z))z = (^y + 37'tan7^ (7' - -0^ sec 7) - y-zxr-^^^ sec 7 (7.112) 
and 

(^5'2)a; = - sin 7 cos-f + 2-f'sm'^-f + {'ip^+w'4)^^)cos-f^smj. (7.113) 
Thus (7.109) is equivalent to 

sec 7 - 2j''ajip^^ tan 7 sec 7 - + 67' tan 7^ ip^ sec 7 



+27"(1 -tan2 7) + (7')^(3 - 4tan^ 7) tan7 = , (7.114) 

COS ^ 



which can be written as 



tp^t + 27' tan 7 wip^^ ( 7^ 67' tan 7^ tp^ 

9i"y' 

= 27'(l-tan^7)cos7 + (7')^(3-4tan2 7)sin7 —. (7.115) 

2 cos 7 



Hence 



cos^7^,,„ 2 X , / X 

yj.^ = ; — ^^(-CtJcos 7) + — m(sec7 — tan7) 

7 2 

+ y"[27"(l-tan^7)cos7 + (7')^(3-4tan^7)sin7]di. (7.116) 

for some one-variable function Q'. Therefore, 

cos'^ 0^ 9z^ 

■0 = — Q'(t27cos^7) + — tJ7ln(sec7 — tan7) + (/? 

7 2 

+^ y"[27"(l-tan^7)cos7 + (y)^(3-4tan^7)sin7](it (7.117) 

for some function cp of t. According to (7.108), 

('t'' \ cos^ 'T 9z^ 
— cos 7 + 27' sin 7 ) H 9(ro cos^ 7) H ln(sec 7 — tan 7) + </? 
7' / 7 2 

+ro y"[27"(l-tan2 7)cos7+(7')^(3-3tan2 7)sin7]di, (7.118) 
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COS f 

II = [ --3(ziTcos^7) +ti7 / [27"(l-tan^7)cos7+ (y)^(3-4tan^7)sin7]dt 

7 J 

+— Ci7m(sec7 — tan 7) + (p\ tan 7 + — sm 7 + 7 w, (7.119) 

2 \ 7' cos 7 y 

hx = 3i/ In (sec 7 - tan 7) + [27"(1 -tan^7)cos7+ (7')^(3 - 4tan^ 7) sin7]di 

— + 27' tan 7 tan7H ^^J'(tJ7cos2 7)}. (7.120) 

7/7 

According to (7.91), (7.92), (7.97), (7.103), (7.105), (7.106) and (7.108), 

gt + fj-gz + gg-x -Qi^f 

_ u"y + 2(7')^ \ , 

7'' sin 7 H — —r cos 7 m57 + (7'' cos 7 + 2(7')^ sin + ipt 

\i) ] 

a' \^ a' \ 

+ I — cos 7 + 27' sin 7 I w sec 7+1 — cos 7 + 27' sin 7 1 sec 7 
VY / V7' / 



, 9 J/7' 

+ (7'' cos 7 + 2(7')^ sin7)CT tan7 + j'tp^w H w 

2 cos 7 



6(7" + (7')^ tan 7) sin 7 + ^ ^^y^2'^^ ^ ^os 7 j ^ + + 27' tan 7^ ^ 



9i"y' 

+i/jt + {Y cos-f + 2{-f'y sm^)w + - — —w, (7.121) 

z cos 7 

— sec^ 7 + ^''w sec 7 + 2(7')^tu tan 7 sec 7 + (7')'^z£7 sec 7 

+ + + l^gz - Qi^f) tan 7 + 2yg^ sec^ 7 
= 2[(7''cos7 + 2(7')^sin7)ro + 7''0] sec^ 7 + 7''rosec7 + 2(7')^^tan7 sec7 

+(7')^^ sec 7 + 2[(7'' cos 7 + 2(7')^ sin 7) sin 7 + 7Vro cos 7]^ sec^ 7 

+ {gt + gg-x + /ifl'^ - 6z//) tan 7 
= 3(7" + 2(7')^ tan7)ro sec 7 + ((7')^ + 27" tan 7 + 4(7')^ tan^ -f)w sec 7 

+27' (V' + ^V'^) sec' l+{gt + QQx + A*^z - 6i//) tan7. (7.122) 

Expression (7.115) says 

91/7' ^ 



( ^ + 27' tan 7 ) ■^ + ■04+ (7'' cos 7 + 2(7')'sin7)ti7 + 
y 7 y z cos 7 

(/?' + [7"(3 - 2 tan' 7) cos 7 + (7')' sin 7 (5 - 4 tan' 7)]^ 

-27'ti7V'^ tan 7 - 27V tan 7 + + 47' tan 7^ 99. (7.123) 
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By (7.121) and (7.123), we set 
9 ^ ^ 



/\4 



6 sin 7 (7 + (7 ) tan 7) + 



yv + 2(y) 



cos 7 I sec 7 



^2 



+^[7" (2 tan 7 - 3 cot 7) + (7') ^(4 tan^ 7-5)] 



+27'^'?/^ sec 7 + 



V?' + ( ^ + 47' tan 7 j ip 



Cc7sec7. 



Then 



(7.124) 



(7.125) 



9x^ 9t + IJ'9z + 99x - 6i// 
by (4.31). Moreover, (4.31) yields d^{F{t,w)) = -d^{F{t,w)) cot 7 for any function F of 
t and w. Furthermore, (7.122), (7.124) and (7.126) imply 

A*t + 9fJ'x + fJ'fJ'z - '^h^T - GuK - 
= 3(7" + 2(7')^ tan 7)ro sec 7+ ((7')^ + 27" tan 7 + 4(7')^tan^7)rosec7 

+27'('^ + wipT^) sec^ 7 + (tan7 + cot 7){[7''(3 — 2 tan^ 7) C0S7 

+ (7')^(5 — 4tan^7) sin7]ro — 27'a7'^^tan7 — 27''0tan7} 
= 3(7" + 2(7')^ tan 7)ro sec 7 + ((7')^ + 27" tan 7 + 4(7')^ tan^ 7)^ sec 7 

+ [7''(3 - 2tan^7) csc7 + (7')^(5 - 4tan^7) sec7]ra7 
= 6;s(csc27 + sec^7). (7.126) 

Therefore, 



P 



-p{32;^(csc 27 + sec^ 7) + 



(fi' + ( ^ + 47' tan 7 I ifi 



wsecy 



+27'ro7/'sec7 + — [7"(2tan7 - 3 cot 7) + (7')^(4tan^7 - 5)] 
2 

-7''%' -|- 2(0'')^ 



(7.127) 



modulo the transformation in (2.21). 

By (4.1) and (7.78), we have the following theorem. 

Theorem 7.4. Let 7, ip he arbitrary functions of t and let '^s be an arbitrary one- 
variable function. Denote w — 2; cos 7 — xsin7 and w — 2; sin 7 + a; cos 7. We have the 
following solutions of the three-dimensional classical non-steady boundary layer equations 
(1.3), (1.4) and (1.5): 



u 



— cos 7 + 27' sin 7 ) H '^{uj cos^ 7) + —w ln(sec 7 — tan 7) + 

7' / 7 2 

+w j [27" (1 - tan^7) cos 7 + (7') ^(3 - 3 tan^ 7) sin7]dt - 



97'^ 9 

IT y ' 

4 cos 7 



(7.128) 
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,cos^ 7 

w — 



4 r 

y^^{w cos"^ 'j) + w / [27"(l-tan^7)cos7 + (7')^(3-4tan^7)sin7]dt 

n . XI fl" . ,l + 2sin2 7\ _ 

+— m( sec 7 — tan 7) + (p\ tan 7 + — sm 7 + 7 tu 

2 \ 7' cos 7 J 



^ ,^ 9 97'ti7sin7 , 

-37'^y sec^ 7 - ^ 7.129 

4 COS"' 7 



3 
2 



\'y^ sec 7 — ( — h 37' tan 7 ) y — In (sec 7 — tan 7) 
> VT / 

-^{ y [27"(1 -tan^7) C0S7+ (7')^(3 - 4tan2 7) sin7]o?i 

+ f — + 27' tan 7^ tan7 + ^^^^Sj'(rocos7)}. (7.130) 
V 7 J 7 

andp is given (7.127) via (7.117). 
Finally, we consider 

fy^ + gy^ + fiy + h, 7] = Ty + Cy'\ (7.131) 

where /, g, h, k, t and C are functions of t, x, z with C 7^ 0. First, 

^y^3fy^ + 2gy + K, ^yy^6fy + 2g, ^yyy^Gf, (7.132) 

^yt = Sfty"^ + 2gty + Kt, = Sf^y'^ + 2g^y + k^c, (7.133) 

^yz ^ ^fzy'^ + 2g^y + K,^, ^ fxy^ + Qxy'^ + i^xy + (7.134) 

Vy^T-Cy'"^, Vyy = '^(y~^, Vyyy = -Ky~^, %x = Tx - Cxy~^ (7.135) 

Vyt = Tt - Cty"^ ?72/2 = - Czy~'^: Vz = r^y + Czy~^- (7.136) 

So we have 

CyCyx - iCx + Vz)Cyy + VyCyz = (3/?/^ + 2gy + k) {3fxy^ + 2gxy + Kx) 

-{6fy + 2g){fy + g^y^ + {kx + Tz)y + hx + Gy~') + (r - Cz/"') (3/.?/' + 2g,y + k^) 

= Sffxy' + ^fxgy^ + (3/./^ + 2gg, - 3fKx - QfT, + 3/,r)y2 _ 2{gC + g.Qy'^ 
+2{gxK - 3fhx - gTz + gzT)y + kKx - QfCz - '^ghx + tKz - 3fzC - (f^zy'^, (7.137) 

CyVyx - iCx + Vz)Vyy + VyVyz = (3/y^ + 2gy + k){Tx - Cxy~^) 

-2Cy~Hfxy^ + gxy'' + i^x + Tz)y + hx + Czy~^) + (r - Cy~')(Tz - C.y"') 
= 'ifTxy'^ + '^gTxy + i^Tx + rr^ - 3fCx - 2/xC - '^{.g()xy~^ 

-{kCx + 2Ckx + 3Ctz + TCz)y-^ - 2hxCy-' - CC.Z/"'- (7.138) 
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Hence (4.2) and (4.3) are implied by the following system of partial differential equations: 

fx = {gOz = «2 = 7-^ = {gOx = /io; = 0, = 6i/, (7.139) 

ft + \ggx - ft^x - 2fT, + f,T = 0, (7.140) 
gt + gxK'- grz + gzT ^0, (7.141) 

+ - 6/G - 3/,C + ^Px = Gi^f, (7.142) 

n + TTz - 3/Ca: + ^P. = 0, (7. 143) 

Ct + /^C. + Xi^x + 3Ct, + rC, = 0. (7.144) 

By (7.139), we take 

/i = 0, f^Ht,z), C^6uz + ^lj{t,x), g^aC' (7.145) 

for some two- variable functions 4>,ip and a function a oft. Moreover, (7.144) becomes 

ipt + i^ipx + '2'4'i^x + 12i^^fi;a; + 3'^T2 + ISuzTz + 6i/T = 0. (7.146) 

Modulo the transformation in (2.19) and (2.20), we assume = or ■^a; 0. Note the 
compatibility — Pzx in (7.142) and (7.143) gives 

/C. = Uufz + Cfzz =^ fi^xx - i^fzz = 201//, + &yzfzz. (7.147) 

Furthermore, the last equation in (7.145) says that g x (7.144) + C x (7.141) yields 

a' + 2a(fi;^ + r,) = (7.148) 

Subcase 1. ipx 7^ 0. 

In this subcase, applying d^dz to (7.146) gives 

4:iyKxx + iJxTzz = 0. (7.149) 

So T is a quadratic polynomial in z. The coefficients of z"^ in (7.146) say that Tzz = 0, and 
so Kxx = 0. By the coefficients of z in (7.146), we have 

K = -2p'x, T = P'z (7.150) 

for some function (3 olt modulo the transformations in (2.17)-(2.20). Now (7.146) becomes 

V't - 2/?'xV^^ - /?V = 0. (7.151) 
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Hence 

z^ = e^9(e2^x) (7.152) 

for some one- variable function Q'. According to (7.147), we take / = 0. Furthermore, 
(7.140) gives us g = 0. 

Subcase 2. ip — 0. 

In this subcase, (7.146) becomes 

2zKa; + ?>ZT^ + T = 0. (7.153) 

So 

= -2/3', T = 7^-^/2 + 13' z (7.154) 
for some functions (3 and 7. Moreover, (7.148) yields 

a = a5^,oe^^ a e M. (7.155) 

Now (7.140) becomes 
and (7.147) gives 

10/, + 3z/,, = 0. (7.157) 
These two equations force us to take 7 = and 

/ = + c, 6,ceR. (7.158) 
By (4.1) and (7.131), we have the following theorem: 

Theorem 7.5. Let (5 be arbitrary function of t and let a, b, c be any real constants. 
Suppose that S> is any one-variable function. We have the following solutions of the three- 
dimensional classical non-steady boundary layer equations (1.3), (1-4) (ind (1.5): (1) 

u^-2(3'x, V ^ (3'y-6uy-\ w ^ (3'z - [6uz + e^Q{e'^^x)]y-'^, (7.159) 



p = p{p" - 2{pr)x' - ^{P" + {pr)z'. (7.160) 



(2) 



2/3 

u^3{be'''^/h-^'^^ + c)y^ + ^^-2f3'x, v = p'y-6uy-\ w = {f3' - 6uy-^)z, (7.161) 

Ol/Z 

p = p[A2ucx + {(5" - 2{l3'f)x'\ - |(/?" + {P'f)z\ (7-162) 
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